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ABSTRACT 

An explanation is developed for the formation, near midnight at 
midlatitudes, of a broad electron-density layer extending approximately 
from 120 to 180 km and usually referred to as the intermediate E layer. 

The responsible mechanism is believed to be converging vertical ion drifts 
resulting from winds of the (2, 4) - mode of the solar semidiurnal tide. 
Numerical solutions of the continuity equation appropriate to the inter- 
mediate layer is described for particular models of ion drift, diffusion 
coefficients and ionization production. Analysis of rocket observations 
of the layers show that the ionization rate is highly correlated with the 
planetary geomagnetic index, K^. Particle flux measurements support the 
idea that energetic electrons are the principal source of this ionization. 
A semiconductor spectrometer experiment for investigation of the particle 
flux, spectrum, and angular propexties was designed and successfully flown 
on a Nike Apache rocket. A detailed description is made of the theory, 
design and calibration of this experiment and some preliminary results 
presented . 
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1. INTRODUCTION 

This report describes in detail the development of rocket-borne 
instrumentation for investigation cf the flux, spectrum, and relative 
directional properties of energetic electrons in the ionosphere. The 
need for these measurements derives from the lack of particle informa- 
tion at midlatitudes under various ionospheric conditions. The signifi- 
cance of particle precipitation as a source of ionization in the F 
region, particularly the altitude range 120 to 180 km, is explored in 
Chapter 2. 

Present theories generally associate the observed ionization to be 

0 

a consequence of ultraviolet radiation (principally Lyman a, 1216 A and 
Lyman 8, 1026 X) scattered from the geocorona and to a minor degree of 
galactic origin. The emphasis here is that energetic electrons may 
play an important role in determining the nighttime E - region iono- 
sphere, particularly at midnight during disturbed conditions. 

In Chapter 3, an explanation is developed for the formation near 
midnight of a broad electron -density layer extending approximately 
from 120 to 180 kilometers. This layer is referred to as the inter- 
mediate layer; it is rbove the region of sporadic-F layers and below 
the F region. In the period from midnight to sunrise the layer is 
observed to gradually narrow and descend at a rate of about 1.5 km min” 1 
until it merges into the sporadic-F layers. The associated mechanism 
responsible for this behavior is believed to be converging ion winds 
(from horizontal wind shears) acted upon by the tidal winds of the 
(2, 4) - mode of the solar semidiurnal tides. The layer is also more 
noticeable during quiet geomagnetic conditions than disturbed. 
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The main part of this report concerns the development of a rocket- 
b 'me solid-state particle spectrometer. The design was developed for 
use on a Nike Apache rocket using IRIG telemetry. More elaborate pro- 
cedures for data manipulation and reduction are available options to the 
system as familiarity, feasibility, and scientific objectives dictate. 
Chapter 4 is devoted to an explanation of semiconductor detector opera- 
tion, sensitivity to internal and external conditions, and comparison 
with other particle detection schemes. Chapter S describes the design 
and fabrication of the spectrometer. 

Discussion is given to optimal design techniques and methods and 
models to be used for interpretation of results. The calibration pro- 
cedure, invoking various experimental procedures, is explained in 
Chapter 6. 
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2. ENERGETIC ELECTRONS AS AN IONIZATION SOURCE IN THE 
NIGHTTIME E REGION 


2.1 Ionisation Rates Calculated from Intermediate Layer Eleatron-Density 

Profiles 

The intermediate layer is observed near midnight in the altitude 
region from 120 to 180 km with peak near 150 km. It is formed as a 
result of the neutral winds, as developed in th? next chapter. Applying 
the continuity equation described in Section . . 5 and integrating 

between the altitude limits z ^ and z ^ gives 



{q - a N*)dz - [jl/zJJ 



0 


( 2 . 1 ) 


where q is the production rate, a is the recombination coefficient, D is 
the ambipolar diffusion coefficient and w is the ion wind velocity 
function. 

Now, if z^ and z ^ are taken to be the nodes of che electron- density 
profile, then dN/dz and w = 0 at z^ and z 7 requiring the square 
brackets in equation (2.1) to be zero. Therefore 


ft 


q dz = 


clN^ dz 


( 2 . 2 ) 


Thus, the height-integrated ionization rate in the layer may be 
obtained without assuming any model for the vertical ion drift velocity. 
The proper choice of limits of integration has enabled the vertical ion 
drift (whether caused by neutral winds and/or electric fields) and 
diffusion to be averaged out; they have not been neglected. 
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An average ionization rate, q, is immediately obtained us 




[aN 2 )dz 


(2.3) 


The integration limits are met at both lower and upper valleys ana 
at the peak of the layer. Thus the average ionization rate may be com- 
puted for (1) the layer below the peak, (2) the layer above the peak, 
and (3) the whole layer. 

The choice of the re coni) in at ion coef cient, a, is dealt with in 
Section 3.2.2. The formulation obtained (equation 3.1) is a function 
of altitude and magnetic index. 

There are now five observations, shown in Figure 2.1, of the inter- 
mediate layer at Wallops Island near midnight: April 12, 1963 

[Cartwright, 196-]; Jane 22, 1965 and February 22, 1968 [Smith, 1970]; 
and September 11, 1969 and November 1, 1972 [Smitn et al . , 1974]. 

Applying the above procedure, the average ionization rate, q, is 
obtained for each of the five cases: for the lower and upper parts of 
the layer, and for the whole layer. These values are given in Table 2.1 
which also includes the altitudes of the valleys and the peaks used as 
limits in the integration of the continuity equation, and the associated 
values of K^. The similarity of the ionization rates over the upper 
and lower portions of the layer suggest that the production is nearly 
uniform over this region. 

The ionization rates are also strongly correlated with magnetic 
index, K^, as shown in Figure 2.2. A least-squares fit to the data 
gives the experimental relationship as follows: 


i 


r 


j 


i 


) 


i 


5 




2 7 


o g 



(««o aonxinv 




u. 


eu 


o 


r— 4 


r-l 




S 


+-» 


c5 




rC 


C* 


•H 


C 


-a 


*H 


e 






<D 


C 


c 


<D 


PS 


cd 

• 


tn 


C 

m 

o 

p 

•H 


+-> 

Q> 

•H 

6 "O 

(1) 

C 


o 


u 

LO 


ai 

u 

0 

•H 

6 

p> 


CD 

<U 

C 


on 

•H 


4«i 

6 

O 

o 


CD 

Cu &o 

X to 

+j 

3 

•H 

o 

to 

• H 

c 

*-» 

<d 

rt 

-3 

> 

i 

u 

o 

CD 

M 

T3 

P 

C 

o 

P 

<u 


*—4 

* 

CD 

as 


•H 

tH 

C 

o 

•H 


e»c 

p 

Jh 

o 

• p 

•H 

> 

P 


rt 

* 

i—4 

"3 


c 


aS 

a 

r—4 

o 

LO 

u 

•— 4 


CM 

0) 

f-i 

3 

oc 

■H 

OU 



Ionization Rate in the Upper E Region near Midnight. 
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I Figure 2.2 Values of and q + , electron production rates in the lower 

and upper parts of the intermediate layer, respectively, as 
inferred by applying equation 2.3 to electron-density profiles 
2.1a-2.1e, as a function of Kp. Smaller dots with indicate 
smaller dots with indicate q * , and larger dots indicate 
q over the entire layer. 
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q = 0.14 exp(0.71 K^cm”'* f*c~* 


(2.4) 


The correlation coefficient for the five observations is 0.96. In- 
terestingly, the data indicate that the lower profile ioniza'icn rate, 
Q^, is always larger for K p values less than 3 while the upper is lways 
greater for K p greater than 3. The least-squares fit corresponding to 
Q l > and Q v are 


Q l = 0.176 exp (0.64 K ) 


Q v = 0.124 exp (0.766 K p ) 


(2.5) 


This reversal suggests that the particle spectrum and pitch angle distri- 
bution is such as to ionize more of the upper portion of the layer as K p 

increases. The greater ionization rate observed at low K for the lower 

P 

profile is expected since Ly 6 ionizes to a greater extent the lower 
E region. Similar results have been obtained by Gough and Collin [1973] 
in which they measured the flux of electrons (> 40 keV) at South Uist 
on thirteen rocket flights. Their correlation with magnetic activity by 
a least-squares fit gives n exponent of 0.98 K p and a correlation 
coefficient of 0.79. 

Applying the ionization rate calculation again to five nighttime 
launches on February 22, 1968, a time evolution of the intermediate layer 
production can be investigated. These results are shown in Figure 2.3 
for the electron- density profiles of Figure 3.1. From midnight to at 
least 0130 the upper and lower ionization rates are similar and quite 
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Figure 2.3 Average ionization rates calculated from the intermediate layer according to equation 
2.3 for the nighttime series of profiles shown in Figure 3.1. The solid dots indicate 
the bottom profile production rate and the plus sign for the upper. 
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small. At 0300 and thereafter there is a substantial increase in upper 
ionization rate presumably from the increasing Ly g flu* from the geo- 
corona and the descending of the layer to lower altitudes. At sunrise 
the ionization quite naturally responds rapidly. 

2.2 Observation of Particle Precipitation in Midlatitudes 

Experimental observations are essential if energetic particles are 
to be recognized as a possible midlatitude ionization source. To 
ascertain the existence of particles, various direct or indirect methods 
can be applied. Particle detectors onboard satellites or rockets pro- 
vide the best direct access to particle information. However, they are 
limited to mappings over a small time scale and require extensive faci- 
lities if simultaneous geographic measurements are desired. Ground-based 
observations rely on measurement of particle aftereffects invoking such 
candidates as light emission, radio noise, magnetic field fluctuations, 
and related ionospheric interactions. 

2.2.1 Rocket-borne particle observations. The most direct measure- 
ments of energetic particles can be made utilizing a rocket-borne 
experiment. This type of experiment is capable of providing excellent 
height resolution, spectrum, and pitch angle information. Although many 
rocket measurements have been made at high latitudes relatively few 
have been concerned with lower latitude observations. In fact, much of 
the present low latitude data is derivable from the results of rocket- 
borne X-ray experiments which are additionally sensitive to energetic 
electrons. 

Two particle detector experiments of the Aeronomy Laboratory were 
launched from Wallops Island (38°N, 75°W; geomagnetic latitude 49°N) at 
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0003 EST on November 1, 1972, and at 2330 EST on April 18, 1974. The 

first payload was equipped with a Geiger counter sensitive to electrons 

greater than 70 keV while the second was equipped with the solid-state 

detector experiment described in Chapters 4, 5, and 6. Both payloads 

included tip mounted Langmuir probes [Smith, 1969] and a propagation 

experiment [Meohtly and Smith, 1970] for electron-density measurements. 

The Geiger counter operation is discussed in Section 4.2.4. The 

particular Geiger counter employed used a beryllium window of area 
2 -3 

0.203 cm and thickness 5.1 x 10 cm. The geometrical factor is 
2 

0.303 cm ster. The window of the counter is oriented perpendicular to 
the spin axis of the rocket and is covered, during the launch phase, by 
a door carrying a weak radioactive source. The door and source are 
ejected 41 sec after launch (43 km altitude) . 

The count rate measured during the flight is shown in Figure 2.4. 

Each data point is the average for a five-second interval. The uncer- 
tainty in count rate is indicated by error bars at representative points. 
The data below 60 km on the descending trajectory are of poor quality 
and are not included in the figure. 

The count rate of the detector at sea-level (before launch) is 
3.3 sec *, of which the source on the door contributes 3.0 sec * and the 
cosmic-ray background 0.3 sec”*. The dashed line in Figure 2.4 connects 
points which have been corrected for the radioactive source carried in 
the early part of the flight. The resulting variation with altitude up 
to about 80 km, including the Pfotzer-Regener maximum near 20 km, is 
consistent with previous observations of cosmic-ray fluxes at mid- 
latitudes [Van Allen, 1952]. In other (daytime) flights with identical 
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counters the background count rate, due to cosmic, rays, has been found to 
be about 10 sec"* at altitudes greater than 40 km. 

The excess counts on this flight are attributed to electrons with 

energies greater than 70 keV. The flux increases with increasing altitude 

in the range 80 to 140 km both on ascent and descent. The count rate 

reaches a maximum value of 339 ± 8 sec"*, which, when corrected for 

-2 -1 -1 

cosmic rays, gives an average flux (> 70 keV) of 1086 ± 26cm sec ster . 

Two important considerations relative to the count rate profiles 
are (1) the detector is predominately looking at particles with pitch 
angles near 90 degrees, but, is modulated because of the rocket spin 
(7.3 rps) and processional (20 sec) periods, (2) the launch was made on 
a very disturbed night, with = 8. The first consideration explains, 
to a large degree, the apparent count rate fluctuation since it is 
observed that the number of particles per incremental pitch angle 
(angular probability) varies appreciably near 90 degrees. The second 
consideration predominately governs the absolute magnitude of the count 
rate curve. 


The launch on April 18, 1974, was made during a moderately disturbed 
night of Kp = 5+. Five channels were used for spectrum analysis corres- 
ponding to energies of 12, 25, 40, 60, and 80 keV. A preliminary study 
indicates that the spectrum near 160 km is related by a power law of 


exponent y - 2 and that the pitch angle distribution varies markedly 


near pitch angles of 90 degrees. Figure 2.5 shows a portion of the flight 
chart record illustrating the typical experiment behavior. The spectrum 


is shown as the staircase waveforms (16 counts per pulse) with the high 


count channel (top one) associated with the lowest energy. The pitch 
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Figure 2*5 Four-inch-per-second chart record of particle data received from 
the five channel energetic electron spectrometer (staircase wave- 
forms) . Also shown are the magnetometer signal (each period 
represents one rocket revolution), Langmuir probe and time code. 

The lowest channel (No. 1) is not shown because of the higher count 
rate. The particle modulation is evident in channel two at the 
top of the chart. 
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angle dependence can be verified in the second channel as a modulation 
with the spin period. The sinusoidal signal is from the magnetometer and 
indicates the relative position of the rocket with respect to the magnetic 
field. 

The particle count rate for the first channel [E > 12 keV) is dis- 
played in Figure 2.6. The results clearly illustrate (1) particle flux 
not changing with time; (2) processional modulation corresponding to 
changing pitch angles; and (3) a nearly linear dependence or count 
rate with altitude. 

These data are amplified by three rocket-borne experiments designed 
primarily to look at galactic X-rays, but also sensitive to energetic 
electrons with energies greater than 4 keV. 

Hill et al. [1970] used proportional counters aimed at right angles 

to the spin axis. The collimation pattern was 5° (vertical) and 30° 

(horizontal) FWHM using the standard egg-crate geometry. The counter 

2 

windows were approximately 60 yg/cra Formvar. Two rockets were 
launched from Kauai, Hawaii. The first on May 15, 1968, at 6:41:20 UT 
(at apogee) and the second on May 17, 1969, at 6:56:15 UT (at apogee). 

The L at apogee (158 km) was 1.16. The rocket spin axis was orientated 
in the vertical causing the counters to scan a band of sky centered on 
the horizon. The resultant count rate verses the azimuth angle for 
both flights is shown in Figure 2.7. The maximum rate of change of 
particle intensity occurs at pitch angles of 90° (corresponding to 
azimuth angles of 270° and 90 6 ) supporting the results obtained from the 
Wallops experiments. Further agreement is made with the linear increase 
in count rate verses altitude and is illustrated in Figure 2.8 for the 
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Figure 2.6 Energetic electron spectrometer low energy (12 keV) channel count rate. The 
data clearly show the 1) temporal smoothness, 2) linear count rate versus 
altitude and 3) precessional peaks due to pitch angle variation. Data were 
sampled every second. 
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Figure 2.7 Count rate verses azimuth from an experiment by Hill t al. [1970] 
for two rocket flights made on 15 May 1968 and 17 May 1969. Solid 
circles show counts corresponding to 1-10 keV energy deposited into 
counter and crosses correspond to energy deposited greater than 
10 keV. 
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May 17, 1969, and for a previous flight on November 6, 1968. All flights 

were made with values of K less than 1. 

P 

Hayakatia et at. [1973] employed two proportional counters equipped 

_2 

with polypropylene windows of 100 and 50 yg cm in thickness and effec- 

2 -2 

tive areas of 32 and 2f cm , respectively. The 100 yg cm counter was 

-2 o 

perpendicular to the role axis and the 50 yg cm” counter was ingled 20° 
upward. Each X-ray event in the energy range 0.1 to 3 keV was pulse 
height analyzed. The equivalent electron energy range was deduced to 
ext 2nd from 5 to 10 keV. 

The launch was made at 2100 Japanese Standard Time on September 3, 
1971, from Kagoshima Space Center, 20°N geomagnetic latitude. The L 
value ranged from 1.17 to 1.25 over the rocket altitude range from 200 
to 800 km. The magnetic dip angle was 42° . The count rate verses alti- 
tude and the azimuthal dependence is reproduced in Figure 2.8. The 

count rate, when multiplied by the appropriate angular sensitive factor 

-2 -1 

of the detector, gives the count rate in units of cm sec ster. This 
count rate is comparable to the count rate of Hill for the overlapped 
altitude range. The dashed line above 400 km represents the altitude 
dependence estimated from the trapping lifetime for the ionization loss 
deduced oy Maeda , [1965]. The curve is normalized to the highest 
altitude experimental data point. The lowest energy data differs signi- 
ficantly and is suggested to be the consequence of atmospheric back- 
scatter. The pitch angle data reveal a very different type of distri- 
bution than that observed by Smith or Hill at lower altitudes. This 
anisotropy, however, is in agreement with the theoretical distribution 
(solid curve) of (sin n a) d(cosa) with n = 8. 
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Physically, the electrons are centered around 90 degree pitch angle 
suggesting that the radiation is trapped and is near the mirroring point. 
The smaller the pitch angle the lower the mirroring point and the 
greater the susceptibility to atmospheric interaction. he 200 to 500 km 
pitch angle distributions indicates a nearly isotropic dependence since 
the higher altitude particles with pitch angles near 90 degrees have 
been filtered out reducing the overall count rate. Continuing below 
200 km the distributions by Smith and Hill show a gradual increase and 
narrowing of the distribution for particles with pitch angles less than 
90 degrees. This is attributed to the increased absorption of particles 
as their pitch angle becomes 90 degrees. The time expended per unit 
height interval is increased as the particle approaches the 90 degree 
pitch angle in its downward trajectory, thereby, increasing the likeli- 
hood of an ionization event and particle absorption. 

The third X-ray flight which additionally revealed the presence of 
energetic electrons is reported by Tuohy and Harries [1973], for their 
experiment launched from Woomera, South Australia, (L = 1.75), at 
0030 hours UT on July 10, 1970. The magnetic index was 7-. The 

detector was sensitive to 25 keV electrons and recorded a flux of 5 
-2 -1 -1 -1 

electrons cm sec keV ster when at 90 degree pitch angle. The 
pitch angle distribution is peaked at 90 degrees with a FWHM spread of 
30 degrees . 

The data from the above experiments intuitively agree with ex- 

i 

pected behavior. Sufficient information exists to deduce ionization rates 
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given the pitch a e distribution, energy spectrum and particle flux. 
These results suggest that the observed midlatitude flux of energetic 
electrons are a primary source of nighttime ionization and will be ex- 
plored in more detail in Section 2.5. 

2.2.2 Satellite particle observations . Satellite data compiled 
by Potemra and Zmuda [1970] are shown in Figure 2.9. The average inten 
sities of electrons were made for energies greater than 40 keV and 
during or near a solar-minimum epoch. Although the maximum intensities 
occur at auroral latitudes it is apparent that there is significant 
midlatitude (35° to 55°) flux and overall flux variability for all 
latitudes. 

2.2.3 Indirect measurements of energetic electron precipitation. 
The indirect measurements of precipitating electrons are of unquestion- 
able importance since they not only support the existence of particles 
but verify the degree of interaction with the atmosphere. This inter- 
action results in such processes as ionization, heating, chemical reac- 
tions, and photon excitations. 

Ionization processes culminate in enhanced electron densities. 
Therefore, by monitoring the electron density by either rocket-borne or 
ground-based techniques, it is possible to ascertain the source magni- 
tude and altitude dependence. For example, if particle ionization is 
observed predominately at altitudes of 180 kilometers, then conclusions 
can be made about the energy spectrum as consisting of 1 to 10 keV 
particles and/or the particle pitch angle distribution is maximized 
around 90 degrees. 
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Figure 2.9 Averaged intensities at sllOO km altitude of precipitated electrons 
>40 keV assimilated by Potemra and Zmuda [1970], satellite obser- 
vations are as follows: (1) Injun 1 during June -September 1961 

[O'Brien, 1962]. (2) Injun 3, January 1963 [ O'Brien , 1964], (3) 

Injun 3, February-October 1963 [Fritz, 1967, 1968, and personal 
communication]. (4) Explorer 12, August-September 1961 [O'Brien 
and Laughlin, 1963]. (5) Alouette, October 1962 to January 1963 
[McDiarmid et al. s 1963]. 
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One particularly interesting observation is the relationship of the 
average ionization produced in the upper E region (120 to 180 km) for dif- 
ferent values of magnetic index, K^. This index is a logarithmic scale 
(from 0 thru 9) of the averaged magnetic field fluctuations. Furthermore, 
this Kp variation is directly correlated with particle precipitation 
since the particle trajectory is perturbed by the changing magnetic field. 
Consequently, Figure 2.2 illustrates the high correlation (0.91) of E- 
region ionization and particle precipitation. 

Another technique which has been correlated with electron precipita- 
tion is VLF (3 to 30 kHz) electromagnetic wave reflection off the D 
region [Potemra and Rosenberg, 1973] . This method is based on measuring 
the phase change of a VLF electromagnetic wave as it is transmitted, 
reflected off the E region, and received at a distant point. The effec- 
tive altitude of reflection and thus, the phase is directly related to 
the D-region ionization. The correlation is then made of reflection 
altitude verses magnetograms showing that enhanced particle precipita- 
tion ( E > 40 keV) , causing a magnetogram fluctuation, ionized a greater 
portion of the D region thereby changing the reflection height. 

Airglow measurements in midlatitudes have been strongly correlated 

with magnetic activity [Hirao et at,, 1965] and hence provide an 

excellent particle flux indicator. The particles predominately excite 

the first negative system of N 2 + releasing photons of 3914 X over the 

altitude range 130 to 300 km. Dalgamo [1964] has found the upper limit 

-2 -2 -1 

of the particle flux to be 3 x 10 ergs cm sec assuming that the 
upper bound of the N 2 + intensity is 60 Rayleighs. An excellent take- 
off to this experiment which would provide substantial ionization 


4 

1 

A 


% 


l 



I 


i 


i 


i 


24 

information is to measure the airglow intensity over an extended period 
of time and to relate it with direct particle data from a simultaneous 
rocket-borne experiment. 

The use of balloons to measure X-ray fluxes is another way to obtain 
an indirect particle measurement. Energetic particles colliding with 
neutral constituents create bremsstrahlung X-rays which are able to pene- 
trate to low balloon altitudes (30 to 50 km) . 

2.3 Charged Particle Dynamics in a Magnetic Field 

In order to properly understand tne defining relations and role 
of charged particles as an ionization source, attention must be given to 
the governing equations of motion of a particle in a nonuniform electro- 
magnetic field. Conclusions may then be reached about particle fluxes 
and the energy spectrum; the variation of both with altitude and finally 
the ionization that is produced. A summary will be given of charged 
particle motion in a uniform and variable static magnetic field. The 
study of the variable magnetic field will lead to the fundamental adia- 
batic invariant with particular emphasis on the magnetic dipole field 
approximation to the earth's field and the associated implications. 

2.3.1 Charged particle motion in a uniform magnetic field. The 
general equation of motion for a particle of charge Ze, according to 
Newton's laws, may be stated as the sum of forces due to it, and non- 
electromagnetic forces f , These forces acting on the particle must 
equal the time derivative of momentum ? 

§.Z<r (7x1.1) .F ex (2.6) 
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For the case of a uniform static magnetic field with zero E and F , 

&X 

the energy of the particle must be constant. Hence 


do tr 


( 2 . 7 ) 


where 


B ~ y 

1 M 


( 2 . 8 ) 


and is the relativistic mass and w B is the angular particle velocity. 
By choosing the reference frame such that the unit vector points in 
the direction of the B field the first order matrix equation (2.7) may be 
simply solved, giving 


t = (Vj Ej * 1> 2 £ 2 ) e 


* U 11 e 3 


where v^, and y^ (parallel velocity) are the integration constants. 
Uj and i >2 can be calculated by specifying a circular motion in the 
and ^2 plane. Hence 


y(t) = w B a (c x - it 2 )e 


' “ll e 3 


(2.9) 


where a is the orbit radius. y(£) may be verified to be a solution of 
2.7. For physical meaning the real part of y(t) is chosen. 
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Integration of equation (2.9) yields the trajectory information 


*(*) * X Q ♦ u n tc 3 + ia(e 2 - ie 2 )e 


-w B t 


( 2 . 10 ) 


where X Q is the position vector of the particle at time zero. The path 
mapped by x(t) is a helix of radius a and with pitch angle a, defined to 
be the angle between the field line and the velocity vector of the 
particle. Hence 


tan o = 


"11 

V 


( 2 . 11 ) 


A useful relation is further obtained from equation (2.8) by sub- 
stitution for y uL = "P x r = P a e_. 

m B 3 


-h-** iS 


(2.12) 


where F is called the magnetic rigidity. 

2.3.2 Charged partiale motion in a uniform magnetic and external 
force field. The next case of interest is charged particle motion in a 
combined static external and magnetic force. Equation (2.6) applies with 
~P set equal to Ze$ + f , yielding two component parts: 


~W 



Si 

dt 





ZeU x $ 


where t is in the reference direction. 


(2.13) 
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The first equation is the conventional acceleration of the charged 
particle in the direction. The second equation can further be de- 
composed into two component parts giving 

P x = U 

(2.14) 

P = Fi t - ZeB x 

y 1 z 


The motion executed by the particles in equation (2.14) is a cycloid 

in the plane perpendicular to 5. The motion is further illustrated in 

Figure 2.10 showing a uniform drift, V^, perpendicular to ~P and t 

(sometimes called an drift). The velocity V„ may be obtained by 

t 

taking an average of equation (2.9) resulting in P =0 and 

y 


p _ + F ext 


ZeB 


B 


eZB 


(2.15) 


Positive and negative particles drift in opposite directions for 

P , — 0 and in the same direction for t, finite. For the case in which 
ext -L 

F > ZeB the ^ field is not able to compensate the strong Ft field 

Am 3 JL 

and the particle is accelerated in the F direction. 

If the external force is a function of spatial coordinates, 
equation (2.6) still applies. The drift velocity in this case follows 
the equipotential line since P = -VP where U is the potential energy. 
Thi? is further illustrated in Figure 2.10 for a particle with 90° 
pitch angle in a uniform S field. The ur. Jerlying approximation assumed 
here in the averaging is that the radius of curvature of the equipoten- 
tial line is much larger than the particles gyroradius. 
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Figure 2.10 Perpendicular motion of a charged particle in a uniform 

static magnetic field with an applied force field F. The 
upper figure show^ that^the averaged particle motion is 
perpendicular to B and F. The lower figure illustrates the 
motion in a nonuniform potential force field. 
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2.3.3 Charged partiole motion in a nonuniform magnetic field. The 
spatial gradients in a magnetic field can be divided into component form 
as VB ■ VjjB + Vj|S relative to the magnetic field t. From Maxwell's 
equations, » 0 and "$x? ■ 0, it follows that if VjB i 0 then VjjJ? is 

also nonzero. The ^contribution will be examined first. Since the B 
field is constant in time, no energy is transferred to the particle so 
that u^ remains constant. The gyroradius, a = (mt’j)/(qB) where B is the 
average field enclosed by an oxbit, remains constant. The force exper- 
ienced by the particle in the y direction does not average out over one 
revolution because of the field inhomogeneity and contributes to a y 
directed force as previously discussed. is directed opposite to V^fl 
and is given by the instantaneous y component as 


f (9) * f cos 6 * qv ( B + V.Sa cose) cose 

a -Li 


(2.16) 


where e is the angle of orbit curvature related to the positive y direc- 
tion. To obtain the average y directed force over one cycln, the inte- 
gral is performed over 0 from 0 to 2it. 
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The associated gradient drift velocity, is computed from equa- 
tion (2.15) giving 
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The motion executed is the same as shown in Figure 2.10. Again 

the implicit assumption made is that V^Ba « B or the orbit curvature 

£ 

must be smaller than the field line curvature, =-a- ■ R , 

n O 

Additionally, if the parallel velocity, Vjj, is nonzero (pitch 
angle not 90°) the parallel curvature of the field line contributes to 
a curvature drift where the centrifugal force is directed perpendicular 
to the 3 field line. 
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and the curvature drift is given by equation (2.15) as 
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The vector combination of these drifts V +7 = ^ is 
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where n is a unit vector in the opposite direction of the center of curva 
ture. Reference to the upper diagram in Figure 2.11 depicts the curved B 


behavior. 



Figure 2.11 The upper figure shows the outward drift of a charged particle 
as it traverses in a curved magnetic field. The drift velocity 
is denoted by V c . The lower figure illustrates the principle 
of reflection of a charged part’de in a converging magnetic 
field. This is the reason for why charge particles can be 
trapped in the earth’s magnetic dipole field. 
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2.3.4 Adiabatic invariance in a nonuni form magnetic field . The 

ev 

magnetic moment, u, is defined to be the product of a current Cjjjp) 

2 

enclosing an area it r and hence 


V 


IA 



v r 


where r is given by equation (2.8) so that p becomes 

2 

mv W 
y = TT~ = T 

o o 


(2.2i, 


where V is the kinetic energy of the circular particle motion. 

Now in a non uni form magnetic field the force exerted due to the 
nonuniformity is F = pVB. This force tends to slow down an energetic 
electron as it moves downward in the earth's converging magnetic field. 
If the magnetic field gradient is sufficient, the field line force is 
capable of completely stopping the particle longitudinal motion and 
return the particle back up the field line as shown in Figure 2.11. 
Since the earth's magnetic field lines are convergent on both sides of 
the dipole the particles are trapped in the geomagnetic field by con- 
tinually being reflected back and forth. More precisely, this reflec- 
tion point or mirror point occurs at an altitude Z, where B - so 
that the total change in longitudinal energy is 

j*u 2 ■ \ 1 Mz * f 1 - -c b i - v 

J O 1 0 


and from equation (2.21) 
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the pitch angle a is given by equation (2.11) requiring 
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(2.23) 


If the pitch angle has a magnitude ^ at 2 = 0, where B = B Q , the 
mirror point will occur when 


(2.24) 


sm a. 


The adiabatic invariant can be derived in the following way. Invok- 
ing Faraday's Law, the E.M.F. generated in one orbit of a particle is 


E.M.F. " . || 


The time for one revolution is 2-nr/v and if AB is the associated 
change in B as the particle proceeds downward in this time the time 
derivative is 


requiring 


at = 2 nr 


. 1.2 AB „ AS 
LW =• 2 M b ~ ^ B 


(2.25) 
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} Differentiating equation (2.21) gives 

An = A [ y j = AB = 0 (2.26) 

and is zero by (2.25). 

The adiabatic invariant are obtained by use of equation (2.21) 

p = constant 

/B = constant (2.27) 

2 

Br = constant 

2.4 Particle Spectrum and Pitch Angle Distribution 

The altitude region in which energetic particles ionize the most 
neutrals is determined by the energy spectrum and by the pitch angle distri- 
butions. Figure 2.12 illustrates the sensitivity of the exponent in the 
exponential spectrum. Also shown is the attenuation profiles for mono- 
energetic beams of electrons. 

The power law and exponential spectrums are plotted in Figure 2.12 
for energies from 1 to 10 keV. The energies from 1 to 10 keV are the pre- 
dominant energies for upper B-region ionization and dramatically show the 
sensitivity of the exponent factors E q and X. The spectrums, J , are 
based on the number of particles with energies greater than the energy E. 

The power law spectrum is given by 

J p (> E) = J x E' x (2.28) 
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Figure 2.12 Exponential and power law spectrums are dramatically different for energies less than 
6 keV. These energies are primary for ^-region ionization. The spectrums, are 
given for all particles with energies greater than E . 
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and the exponential spectrum by 

J e i> £) = J q exp i-E/EJ (2.29) 

The E-region particle data have shown power- law spectrums with X 
usually between 1.5 and 2.5. This supports the fact that particle ioni- 
zation is nearly the same in the upper portion of the E region (150 to 
180 km) as it is in the lower (120 to 150) . This is a consequence of the 
ionization of the upper E region by predominately 1 to 4 keV electrons 
and the large increase in 1 to 4 keV particles associated with the power 
law spectrum. 

Consideration must also be given to the fact that the particle spec- 
trum can vary with pitch angle and consequently with the altitude. The 
process of obtaining a spectrum by averaging over a wide pitch angle 
spread must be justified. The incentive for the exponential spectrum 
is that it is Maxwellian in nature suggesting internal equilibrium. A 
plot of the differential number flux divided by the associated energy 
and plotted against energy will give a straight line for a Maxwellian 
spectrum. 

The particle data show the flux to continually increase with alti- 
tude in a linear fashion at midlatitudes for particles with energies of 
at least 10 to 80 keV. The only viable explanation is that equal energy 
particles are being reflected at differing altitudes because of the 
mirroring point variation (equation 2.24). In other words, the pitch 
angle distribution is such as to cause particles of a particular energy 
to have a linear variation in mirroring altitude. Since the parallel 
velocity is zero at the mirror point, equation (2.22) becomes 
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and since the pitch angle is related to the velocities as sina^ 
the previous equation is expressed as 

B 

S(r,0) = j— 

sin o 

o 

Now since an equal number of particles are reflected in each height 
interval over the concerned region, the solution of equation (2.31) for 
will yield the predicted pitch angle distribution for any trajectory 
point invoking the earth's magnetic field variation. 

2.5 Ionization Rates Calculated from Flux and Spectrum Information 

The particle ionization rates may be calculated from knowledge of 
energetic particle flux, spectrum, and pitch angle together with infor- 
mation on the quantum efficiency of neutral constituent interaction. 

Rees [1963], Berger et al. [1970], and Wulff [1973] have investigated 
this problem for cases particularly associated with the higher latitude 
ionosphere. Important differences in particle pitch angle and ionosphere 
composition at night warrant further investigations. 

The techniques developed are modifications of the methods employed 
by the above. The computer program developed for this calculation is 
given in Appendix I . Utilizing the method of Wulff, the ionization 
curves shown in Figure 2.13 were obtained for an isotropic pitch angle 
distribution. The results illustrate the sharp attenuation of particles 
with altitude. The curves are significantly altered by different pitch 
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angle distributions with greater dependence as the energy is decreased. 
This is exanplified by calculations from Prasad and Singh [1972] in 
which they deduced pitch angle dependence (Figure 2.14) of energetic 
particles attenuated in the earth's atmosphere. 








ALTITUDE (km) 


4 



ENERGY (keV) 


Figure 2.14 Pitch angle variation in ionization curves for monoenergetic 
electrons of 4 and 40 keV. Data replotted from Prasad and 
Singh [1972]. 
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3. DYNAMICS INVOLVED IN THE FORMATION OF THE 
E REGION INTERMEDIATE LAYER 
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3.1 Morphology of the Intermediate Layer 

Experimental measurements of electron densities indicate the pre- 
sence of a layer in the upper E region near midnight. The layer usually 
extends from 120 to 180 km and has a peak near 150 km. The time evolu- 
tion of the layer is observed to narrow and descend at approximately 
1.5 km min"* until it eventually appears as a morning sporadic-# layer. 
Figure 3.1 presents five rocket observation [Smith, 1970] made from 
Wallops Island (38°N, 75°W; geomagnetic latitude 49°N) which detail this 
descent. 

Five rocket observations made at midnight (Figure 2.1) of the 
intermediate layer, suggest a strong relationship with K^. The layer, 
at low values of Kp, is quite sharply peaked as compared to the almost 
uniform profiles at high values . 

Geographically, the intermediate layer is observed at middle lati- 
tudes. Extensive studies concerning the global characteristics of the 
intermediate layer are investigated by Wakai and Sawada [1964]. 

The critical frequency of the intermediate layer increases pole- 
wards from latitudes 0° to 50° and decreases thereafter. Additionally, 
the critical frequency of the layer has a positive correlation with 
[Berkner and Seaton, 1949]; [ Watts and Brown, 1954], 

3.2 Continuity Equation as Applied to Vertical Ion Transports 

3.2.1 Hind shear theory as utilized to explain the intermediate 

layer. Explanation for the ionization bunching in the vicinity of 
150 km must reside in some dynamical forcing function since nighttime 
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Figure 3.1 


Electron-density profiles from the probes on rockets launched 
from Wallops Island on 2*. February 1968, [ Smith , 1970]. The 
profiles are each displaced by one decade. The electron- density 
scales for the first and last profiles are given. The part of 
the profiles below 90 km should be regarded as indicating the 
structure, rather than as absolute values. 
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production and recombination are believed uniform. Previous investiga- 
tions suggest that the vertical redistribution of ionization is a conse- 
quence of the neutral winds [Conatantinidea and Bedinger, 1971; 

Fugitaka et al. , 1971]. Basically, the predominant hcrizontal winds flow 
in opposite directions, one above the other, causing a shear in the 
vicinity of ISO km. Now, the motion of ionization is subject not only to 
the force of the wind but also to the magnetic field. The product 
gives a net upward flux and downward flux component dependent on the 
wind direction. For midlatitudes the wind system is such as to cause 
ionization to migrate towards the neutral wind shear region. Ionization 
is therefore transported to the peak from a lower and higher region 
surrounding the peak. An extensive analysis of vertical ion drift 
velocity from neutral winds is formulated by MacLeod [1966]. Electric 
fields may also contribute to the vertical ion drift velocity; the 
following analysis is independent of the mechanism producing the verti- 
cal ion drift, however. 

3.2.2 Production and recombination. The source of ionization has 
been presumed to be solar UV radiation resonantly scattereo by the geo- 
corona, principally Lyman-a (1216A), which ionizes nitric oxide [Suider, 

O 

1965] and Lyman-6 (1026A), which ionizes molecular oxygen [Ogaua and 
Tohmat8u, 1966] . These UV sources are believed to be of secondary 
importance as compared to that of energetic electrons for the reasons 
enunciated in Chapter 2. Also, it has been found [Meier, 1974] that 
radiation is at a minimum near midnight. 

The altitude variation of production for input into the continuity 
equation was initially taken to be constant. This was based on the 
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material presented in Section 2.4. Here ionization rates were deduced 
for the peak to upper node and for the peak to lower node from the inter- 
mediate layer electron-density profile and appropriate recombination 
coefficient. The justification for using a constant production rate 
over the concerned altitude is the result of the close agreement between 
upper and lower ionization rates. Further refinements will take into 
account variable production as data becomes available. The developed 
program for solving the continuity equation is compatible for such 
variations. 

The recombination coefficient is calculated based on the particular 
ions present and their temperature. At night, in the altitude range of 
interest, the principal ion is N0 + [Holmes et at., 1965]. The loss 
process is presumed to be dissociative recombination of N0 + ; the coeffi- 
cient (o) being given by Biondi [1969] to be 

a = 4.5 x 10 -7 (300/T,) cm 3 sec -1 (3.1) 

The electron temperature (T ) above 120 km, however, varies with geo- 
magnetic activity. For purposes of calculation the data from Evans [1973] 
at Millstone Hill are employed such that the nighttime (2100-0300 EST) 
electron tenperature at 200 km can be well represented by (ll2Kp + 586) °K 
for eight occasions in 1966 and 1967 with a range of values of of 1- 
to 5-. The correlation coefficient is 0.90. Following Evans the elec- 
tron tenperature at 120 km is assumed to be T g ■ 355 # K. A linear inter- 
polation for altitudes between 120 and 200 kilometers gives 
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T (a) = (3-120) x (1.40 K + 2.91) + 355°K (3.2) 

e p 


Using Biondi's formula for a and allowing the extrapolation of electron 
temperature to = 8 the recombination rate is obtained as a function 
of altitude and geomagnetic index. 

3.2.3 Solution to the continuity equation and sensitivity analysis. 
The continuity equation for electron density is fundamental to an under- 
standing of ionospheric variations. The basic conservation principle 
gives the temporal change in electron density at each height as 


3 *e 

- q - l • div (t yy 


(3.3) 


where the three right-hand terms represent the ionization production 
rate, the ionization loss rate, and the effect of transport, respectively 
The transport term is made up of movements such as diffusion, electro- 
magnetic drifts, and thermal expansions. 

For the steady-state case with no transport Ratcliffe [1956] has 
formulated the loss term to be 


L = LN 


ill 


N + (-) [I] 

e or 1 J 


(3.4) 


where a and X are, respectively, the rate coefficients of dissociative 

recombination and ion-atom interchange reactions. The brackets denote 

the number density of species X. The associated [X], a, X, and N for 

the E region is such that N « X[J(] for the involved reactions. The 

2 

recombination term is consequently of the form o iN . 
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The divergence term is responsive to the flux motion. Of particular 
interest is the flux motion due to a forcing function giving the elec- 
trons a velocity ^ and the motion resultant of the gradient concentration, 
(diffusion) . The diffusion flux term from elementary deduction is giv»n 
as 


* « fa = D$N 


(3.5) 


The general time dependent continuity equation for electrons in the 
E region, therefore, is given by 

S ♦ $ • (fo) * q ■ atf 2 ♦ $ • Dh (3.6) 


dN 

This equation is simplified by omitting the term (since obser- 
vations suggest the time variance to be small) along with the terms 
arising from the horizontal part of the gradient operator. Permitting 
to be the vertical ion drift velocity, and z the altitude, 
equation (3.6) ^educes to 

= q - otf 2 + j£(Z>ff ) (3.7) 

For the case of interest, the ion wind system as presented in 
Section 3.2.1 is convergent towards the intermediate layer peak. To 
ascertain the behavior of equation (3.7) a sinusoidal wind approxima- 
tion is made. The actual winds measured, in fact, »re in quite good 
agreement with this approximation. Of course, many such sinewave 
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solutions could be solved for and added in a Fourier series for general 
behavior. The wind system, therefore, representative of the intermediate 
E layer is taken to be 


w = 20 sin [it ( 150 -a)/ 30] m s' 


(3.8) 


Here, the amplitude is taken as 20 m sec * with the layer peak at 150 km 
and upper and lower valleys at 120 and 180 kilometers, respectively. 

Appendix II describes the solution of this equation for the case 
D = 0 with the assumed sinusoidal vertical ion drift. r ..e to the 
numerical sensitivity of this equation towards instability a precision 
predictor- corrector Euler-Romberg method was used for analysis. Addi- 
tionally, a complicated analytic series solution was derived for the 
D - 0 case. The no-diffusion continuity equation may be transformed as 
the Riccati equation. Using various substitutions, a second order 
linear differential equation is obtained. Transformation again results 
in a form applicable to the WKB method of solution. 

To solve the variable diffusion continuity equation (Appendix III) 
a Runge-Kutta method was applied with 0.2 km steps. Further, since the 
initial conditions cannot be specified at one point, as required by this 
iterative scheme, a convergence procedure is required. The two known 
boundary conditions of tne second order equation are zero slope (zero 
flux) at the peak and valleys of the profile. Thus, by starting at the 
peak of the profile and assuming a typical electron density (i.e. a 
value between the D = 0 value and the photoequilibrium value) a 
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convergence towards the actual point will occur, if one monitors the 
final slope of the valley and corrects accordingly. 

Solutions with and without diffusion are shown in Figure 3.2 for 
five values of Q. The thick profiles represent a diffusion coefficient, 

g 

D, of 1 x 10 and the dotted profiles for D = 0, the results clearly 
indicate that diffusion is significant for low production ”ates. 

Results illustrate that the modeled nighttime intermediate layer is in 
close agreement with observation (Figure 2.1) and is therefore a plausi- 
ble explanation. 

-3 -1 

Profile 1 (S = 0.4 cm s ) gives results looking quite similar to 
the intermediate layer at midnight under quiet conditions. As the ioni- 
zation rate increases (eg. increases according to Figure 2.2) the 
peak broadens and becomes less pronounced. Figure 3.3 illustrates the 
intermediate layer dependence on using the relation developed in 
equation (2.4). The peak of the layer is similarly related to Kp as 

_7 

illustrated in Figure 3.4 for diffusion coefficient- of 0, 0.5 x 10 
-7 2 

and 10 cm sec. 

The sensitivity of equation (3.7) to vertical ion wind variations 

in amplitude is illustrated in Figure 3.5. The computations are carried 

out for values of amplitude. A, of 0, 10, 20, and 50 m sec * . Much 

ionization redistribution is seen to occur at relative low values of 

the ionization drift velocity. However, this redistribution does not 

affect the electron-density profile in a proportional manner. For 

example, doubling A from 10 to 20 is seen to increase the peak value of 

3-3 3 3 

electron density from about 5.1 x 10 cm to 7.5 x 10 cm , an increase 


of about 50%, while increasing A from 20 to 50 gives an increase of 
only 20%. 
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3.3 Winds Derived from Electron- Density Profiles 

The electron-density profiles of the intermediate E -region layer may 
be used to deduce information of ion and therefore neutral winds . For 
this purpose the continuity equation is employed using the deduced ioni- 
zation rates calculated for the upper and lower part of layer 

(equation 2.3). Referring to the continuity equation (3.7) and taking 

dN 

the lower integration limit z ^ to be the altitude where w = = 0 and 

rearranging, becomes 


w(a) = 


(<7 - aN 2 )dz + D (^) 


(3.9) 


where r is evaluated as in Section 2.1. This expression gives the 
vertical ion wind velocity in a readily available form as a function of 
production, recombination, diffusion, and electron density. The program 
used to calculate the average ionization rate, q, and the vertical ion 
wind is presented in Appendix IV. The program requires electron den- 
sities every 0.5 kilometers. The recombination (Section 3.2.2) and 
diffusion coefficient are taken to be variable. 

3.3.1 Winds calculated for various disturbed conditions. Solving 
equation (3.9) for the intermediate layers (Figure 2.1) measured during 
Kp values of 0+, 1, 2+, and 3+ yields the associated wind patterns of 
Figure 3.6. The sinusoidal wind pattern is apparent. No strong correla- 
tions exist; however, some justification can be made for a positive 
correlation of the wind root-mean-squared-amplitude verses K^, No 
reasonable hypothesis for this dependence is evident. 
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3.6 Vertical ion wind velocity calculated from electron- density profiles for 
values of 0+, 1, 2+ and 3+. 
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The lower amplitude deduced for the upper ion wind has not been 
fully evaluated. Assuming the approximations made are indicative of this 
region, the straightforward conclusion is made that the velocity of the 
upper wind system is slower than the lower. A high correlation exists 
between the ratio of the upper and lower amplitudes. 

The profile for = 8 is not shown because of the large amplitude 
variation. During high values of K^, the ionization rate is not believed 
to be uniform as for lower values. Consequently, erratic behavior is 
apparent in the electron-density profile for this case. The further 
incentive for not including large wind amplitudes is due to their inhe- 
rent lower sensitivity to electron- density profiles as indicated in 
Section 3.2.3. 

3.3.2 Winds calculated from a nighttime series of profiles. On 
February 22, 1968, five rockets were launched at 1.5 hour intervals 
beginning at 0009 EST. The electron-density resrlts (Figure 3.1) 
clearly show the intermediate layer descending at approximately 1.5 km 
per minute. The winds derived from these profiles are shown in 
Figure 3.7. The assymetry of the ion wind profiles changes immediately 
as the layer begins to descend suggesting that the relative wind 
velocity reverses. 

The descending layer can be associated with a tidal wind fluctuation. 
The apparent vertical wavelength (about 50 km) and nodel descent rate 
imply a wave period of about 14 hours. This is not too different from 
the parameters of the (2,4) -mode of the solar semidiurnal tides, the ver- 


tical wavelength being around 50 km and the period being 12 hours. This 
conclusion has likewise been reached by Fujitaka and Tohmatsu [1973], A 
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theoretical prediction of an enhanced (2,4) -solar semidiurnal tide above 

100 km has been put forth by Lindzen and Hong [1974] . 

The erratic behavior of the » 8 electron- density profile and the 

0602 early morning profile are used to derive the wind system shown in 

Figures 3.8 and 3.9. Both profiles are substantially outside of the 

limits of the original assumptions requiring that further justification 

be made. The erratic behavior of the K = 8 profile leads to the erratic 

P 

derived wind system of Figure 3.8 if a dynamical explanation of wind 
shears is assumed. Irregular ionization is also a possible candidate 
due to erratic spectrum and pitch angle variations. The low (0602) 
profile needs to be re-evaluated for the appropriate recombinations 
coefficients associated with the lower E region. The rapid shear 
illustrated in the figure suggests an unnatural wind system and therefore 
additional forcing functions (i.e. electric fields) must be likewise 
investigated. 
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VERTICAL ION DRIFT VELOCITY (msec” 1 ) 

Figure 3.8 Anomalous vertical ion wind velocity calculated from a K = 8 electron -density profile. 







Figure 3.9 Anomalous vertical ion wine 1 velocity calculated from an early morning (0602 EST) measurement 
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4. THEORY OF SOLID-STATE PARTICLE DETECTORS 

4,1 Quantituvive Theoretical Consideration 

Energetic charged particles propagating within a solid lose their 
kinetic energy through lattice interactions. These interactions may be 

■rtf t 

approximately thought of as similar to ionization of a gas in a Geiger 
counter except that the gas is in a condensed form, a solid. In a gas, 
the decaying energetic particle r roduces electron-ion pairs which are 
subsequently swept out by an electric f'.eld. Accordingly, they 
increase their energy, yield multiplying secondaries, and together con- 
tribute to a current pulse. 

The correct treatment for a crystal lattice, since there are inter- 
action fields, is to formulate a many-body wave-function of the 
Hartree-Fock type, and then to apply an energetic particle wave as a 
perturbation in the Quantum Mechanical description. For interactions 
which are associated with particle detection it suffices that the "con- 
densed gas" model applies utilizing Fermi-Dirac statistics. Here, the 
allowed energy states (called an energy band for a series of adjacent 
states) depend on the nature and spatial configuration cf the molecules. 

The probability that a particular energy state (level) will be filled is 
given by the Fermi equation and is a function of energy, temperature, and 
impurity concentration. For a semiconductor crystal the energy states 
(from the Hartree-Fock equation) fall into two primary bands (valence 
and conduction) separated by a forbidden energy band gap. Energetic 
particles, impinging on the detector, supply energy to the lattice elec- 
trons lifting them from the nonmobile filled valence bands into the 
essentially empty conduction band. 
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Applying a gradient potential within the material, an electron-hole 
current is established. The hole current is the name for the process of 
sequential exchange of valence electrons between adjacent lattice sites. 

The energy lost in ion-electron pair formation in a gas is approxi- 
mately 35 eV whereas for silicon and germanium hole-electron formation 
is 3.5 eV and 2.9 eV, respectively. The actual band gaps for silicon 
and germanium are 1.1 eV and 0.67 eV. The differential energy discre- 
pancy between the band gap and average pair generation energy is recon- 
ciled by coupling of electrons into lattice vibrations. This is analo- 
gous to the pair production within a gas where surplus energy is taken up 
in excitation and dissociation processes. 

Varii. is other crystalline solids have been used successfully. 

Diamond and cadmium sulfide have bee • used at room temperatures while 
halides of silver and thallium are suitable at low temperatures where 
their ionic conductivity is minimal. Typical charge separation energy 
for these solid dielectrics is 10 eV. Material impurities contribute to 
nonuniform collection of electron-holes pairs and, consequently, limit 
the use of solid dielectrics. This inhomogeneity is called "trapping" 
since the electrons or holes may be attached to high charge (polariza- 
tion) centers. Semiconductors are less sensitive to "trapping" because 
electrons and holes have longer lifetimes and greater mobilities. 

4.2 Comparison of Particle Detectors 

4.2.1 Advantages and disadvantages of a solid-state detector. The 
principle advantages of solid-state particle detectors over gaseous 
devices are summarized as follows: 
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(a) Better statistics and resolution result since many more charge 
carriers are released for a particular incident energetic 
particle. 

(b) Very short deadtimes on the order of a few nanoseconds are 
possible because of the small collection distances and large 
carrier mobilities. The collected charge is independent of 
the location of the ionization event. 

(c) Increased experimental certainty is possible. Problems 
associated with the gas purity and time stability are eli- 
minated as negligible contributions. 

(d) Detectors are rugged, compact, and suited for easy mounting. 

(e) Lower bias (about 100 volts); requires no elaborate corona- 
breakdown shielding precautions. 

(f) Better control is possible of detector sensitive depth, area 
and geometry. 

(g) Lower energy particles can be detected. Furthermore, the 
separation charge energy is independent of the incident 
particle type (electrons, protons, a particles, and heavier 
ions) which results in a linear function of initial energy. 

(h) High stopping efficiencies of particles within a small region 
are possible. The total absorption peak efficiency is 
approximately 0.7 times the geometrical efficiency. 

A solid-state detector suitable for rocket payload applications is 
shown in Figure 4.1. Two primary problems with solid-state detectors 
arv channeling of energetic particles between crystal planes, and the 
relatively weak signals requiring extremely sensitive electronics. 
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Figure 4,1 Ortec rugged i zed surface barrier detector. The 

sensitive area is 50 mm-, depletion depth of 100 um 
operating at a reversed bias of 120 volts and a 
noise width of 6 keV at 25°C. The flat, mirror, 
surface on the right is the sensitive material. 
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These factors are discussed in later sections. The energy range of the 
solid-state detector is compared in Figure 4.2 with other detection 
schemes. Some characteristics of the other principal rocket-borne parti- 
cle detectors are summarized below. 

4.2.2 Channeltron particle detectors. The channel multiplier is 
the most attractive instrument for overall low energy coverage as indi- 
cated in Figure 4.2. The geometrical construction of the linear and 
helix channel multiplier are displayed in Figure 4.3. Figure 4. 4. shows 
a curved channeltron interfaced to an electric-field scanning energy 
analyzer (black section) and collimator. Operationally, energetic 
particles impinging upon the inside surface of the channeltron tube 
opening collide with the resistive wall creating secondaries which are 
subsequently made to accelerate axially along the tube by an applied 
electric field. Further multiplication arises as the accelerated secon- 
daries collide with the walls initiating tertiaries, etc. Typical gains 
9 8 

are of 10 in linear and 10 in curved channeltrons. 

Complications in the application of channel multipliers are 1) the 
gain is pressure sensitive due to the ionization of the neutral residual 
gas (ionic feedback), 2) necessity to vacuum seal the tube to prevent 
outgassing, 3) high voltages required to induce the electric field, 
and 4) early saturation begins at 10^ pps because of field distortion 
and space charge effects. The chief reason for implementing a curved or 
spiral channeltron is to reduce ionic feedback. Further reference can 
be made to the following publications by Adams and Manley [1966], 

Reed et al. [1969], and Sacerlotzky and Belanger [1972]. 
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Hie linear channel tron in the top photo was used for cali- 
bration of the solid-state detector discussed Ln Chapter (>. 
Particles passing through the grid collide with tin walls 
and experience an applied axial electric field causing far- 
ther collisions and multiplication. The bottom > into illus 
trstes a rocket-borne helical channeltrun used to reduce 
pressure sensitivity and overall dimension. The horizontal 
lines are 0.8o cm apart. 
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f'i t ,urc 4.4 The curved channel t run, coated in epoxy, is attached to an 

electric field scanning electrode used for spectrum determina- 
tion. The entrance and collimator section precede the electric 
field analyzer and provide the proper particle angular pattern. 
Properly calibrated and used in conjunction with a so.i i- state 
detector, complete, accurate and high resolution particle cover- 
age is obtained. 
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4 . 2.3 Scintillator particle detectors. Energetic particles 
colliding into a suitable solid or liquid material decay exciting elec- 
trons which in turn release photons (scintillations). The light is 
detected in a photomultiplier which is calibrated to give the proper 
particle energy. Primary advantages of the scintillator are its sim- 
plicity, large sensitive volume, and complete spectrum readout capa- 
bility in a counting interval unlike that of a sweep magnetic spectro- 
meter. The scintillator is applicable to energies down to approximately 
2 keV; however, the main drawback is that there is poor energy resolu- 
tion. The scintillator response to monoenergetic electrons is essen- 
tially a gaussian peak whose width varies inversely with the square 

root of photocathode and multiplication electrons in the photomultiplier. 

4.2.4 Gaseous ionization detectors. The gas detector is classi- 
fied into three modes of operation: ionization, proportional, and 

Geiger operation. Figure 4.5 is a closeup view through the window of a 
Geiger counter used to detect X-rays and energetic electrons. Counting 
data are shown in Figure 2.4, for this unit. The collection scheme is 
based on the principle of the ionization of a gas by a charged particle 
whereby the ions are accelerated in an electric field yielding a current, 
pulse. If no secondaries are formed (low electric field) the output 
pulse is simply the very small charge released by the decaying particle 
and is proportional to its initial energy. When used in an integrating 
mede this is an ionization chamber. The proportional counter utilizes 
the multiplication of secondaries by applying a higher electric field 
strength. Increasing the field still further gives complete breakdown 
of the gas in the chamber and is termed the "Geiger region". Although 




Figure 4,5 Displayed is a proportional counter chamber. 

Energetic particles pass through the window to 
ionite 1 he enclosed gas. The applied electric 
field accelerates the ion iced charge causing 
further ioni ratio is and increased current flow. 
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a simple experiment, the gas detectors suffer because of 1) gas purity, 
2) position sensitiveness with respect to where the ionization track 
takes place within the gas, and 3) associated deadtimes due to slow ion 
mobilities. 

Detection schemes which have been employed to a minor extent are 
photographic film, spark chambers, spacecraft charging [DeForest, 1972] 
and tracks left in plastic, glass, or micasheets. A review of the 
various detection methods may be found by O'Kelley [1962]. 

4.3 Diode Geometry and Construction of Solid" State Detectors 

Signals introduced by an energetic particle are usually very small 
CIO -15 coulomb) requiring every precaution to insure maximum signal -to- 
noise ratio. The reverse current flow with the applied electric field 
appears as noise at the output. It is minimized by increasing the 
effective resistivity. One of the most attractive procedures to 
accomplish this is the p-n semiconductor diode junction in which large 
electric fields can coexist with low reverse leakage currents. A 
qualitative analysis will be given in Section 4.4 of diode behavior. 
Other techniques used to increase resistivity are the implantation of 
gold impurities at low temperatures [Davis, 1958] and high temperature 
imperfections introduced by temporarily heating crystals to 1000°C. Two 
principal problems associated with these methods are induced "trapping^ 
center: and necer-sity for low temperature hardware, which mo’ e their, 
unsuitable for space applications. 

Figure 4.6 illustrates the basic geometry of three types of semi- 
conductor detectors; the diffused (type I), surface barrier (type .1), 
and the lithium drift (type I T I). The type I . 'nsitive surface consists 
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Figure 4.o Three types of solid-state particle detectors are represented showing their respective 
geometry and electric field strengths. For totally depleted detectors, the depletion 
region depth is equal to the total detector length L. 
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of a thin highly doped n+ donor impurity while the type II case is a 
metallic surface e.g. gold in contact with the n material. Two quali- 
fications placed on the diode are that 1) it must be fabricated such that 
the depletion region extends over a range somewhat larger than the inci- 
dent particle decay path if energy spectrum information is desired and; 
j 

l 

2) energy lost in the thin surface region (dead zone) before entering 
I the depletion region is optimized to its smallest reasonable valve since 

. t does not contribute to the output signal. In the depletion region 
an electric field exists resulting from the redistribution (diffusive 
process) between adjacent materials of the unequal number density of 
conducting band electrons and free holes (maximum entropy occupancy of 
energy and space states). 

For large energies requiring wide depletion regions the Lithium ion 
I drift (type III) technique was introduced by Pell [I960]. Layers 

approximately 15 mm long are possible with nearly intrinsic resistivity 
at low bias voltages of 100 V. Firstly, this method, like normal junc- 
tion fabrication, allows diffusion of donor (lithium) atoms into the 
crystal by passing a high concentration of donor atoms at elevated tem- 
peratures over the surface of a low-resistive p-type crystal. A junc- 
tion is formed at the point where the concentration of ^ acceptors 
equals -hat of n donors. Secondly, to spread the depletion width Pell 
applied a reverse bias to the junction at conti iued elevated tempera- 
tures causing lithium ions to jrift into the p-type base under the action 
of the electric field. The charge-compensating effect * f the Lithium 
donor movement in the electric field offsets the acceptor impurities to 
give a uniform region, the i layer, of intrinsic resistive mate" 1. 
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4.4 Qualitative Considerations of Particle Impingement in Type I and II 

Diodes 

To properly interpret the retrieved data from a solid-state detector 
it is of fundamental importance that discussion be made of the basic theory 
of operation and the sensitivity of the detector to external conditions. 

This section is intended to develop an intermediate working language of 
the detector operation with emphasis on the pertinent factors involved in 
the rocket-borne design of Chapter 5. Advanced and detailed treatments are 
in the publications by Deme [1971], Bertolini and Coche [1968], and 
Taylor [1963]. 

4.4.1 Energy deposition in solid-state detectors. Figure 4.7 shows 
dE 

relationships for electrons replotted from Stermheimer [1959] and for 
protons and alpha from [Williamson et al., 1962]. Channeling of ions between 
crystal planes can cause deviations from the data. The energy loss is 
primarily due to Coulomb interaction of the incident charged particle with 
excitations of the absorber electrons. The functional dependence has been 
deduced by Livingston and Bethe [1937] to be 


dE 

& = ’ 


4 2 

47T<? z ZN 



In 



[ erg/cm ] 


(4.1) 


where E is tl e kinetic energy of the incident ion, x is the particle track, 
e is the unit charge, z is the charge of the particle, and 2 is the atomic 
number of the absorber material, N is the absorber density (cm , m is 
the electron mass, v is the particle velocity and I is the average 
ionization potential of the absorber. 
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Letting M be the mass of the particle, C4.11 becomes 


% “ T \ ln y2 + Qomt ] 

V 


z 2 M 


In tt + oorcst 
M 


(4 2) 


Therefore 


E « Ms 2 
dx 


\ in i 


+ const 


- Mz (const) 


(4.3) 


Accordingly, equation (4.3) shows that particle identification is 

dE 

possible if simultaneous measurements are made of (Figure 4.7) and 

total energy E. Iheir product may be formed in a multiplying circuit 

2 

yielding discrimination of z m. A working arrangement employed is to have 
to detectors, a thin total ’y depleted detector of length dx directly above 
an absorbing detector for final energy measurement. 

The range in silicon, R, is found by integrating the specific 
energy loss, equation (4.1) 



dx 
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where E and v are the initial energy and velocity of the particle. The 

o o 

range in silicon for electrons, protons, and alpha particles is plotted 
in Figure 4*8 from data by Williamson et al . [1966], 

4.4.2 Charge signal generated tn a o lid- state detector . The 
shape of the electric field depends on the diode type and compos .ti^n. 


i 



Specific Energy Losses- (Key-micron*' of Si) 






PARTICLE ENERGY IN Kev 





1 


} 

i 


1 


i 


Using diffused and surface barrier detectors, the electric fields are 
nearly linear as shown in Figure 4.6. 


2 V T 

e = ± =-j (L-x) + 

L 


~r u 


where V T is the voltage required to totally deplete (L=<i) the detector. 

dF 

Given information on the electric field (equation 4.5), ^energy 
loss (Figure 4.7) and mobilities, y(e,t) (Figure 4.9) it is possible to 
calculate the charge transit time which together with the "plasma time" 
specify the output charge signal Q(t) . 

To develop this relationship necessitates the use of a fundamental 
theorem first enunciated by Ramo [1939]. The theorem states that the 
induced charge A Q derived from a carrier charge q' transit in a region 
LX (normal component) between two parallel plane electrodes separated 
by a distance L is: 


A n , A* 

A« = — 


Now for charges generated along a track of length R, the above 
theorem suggests that the charge produced in each segment kX can be 
summed to yield the total induced charge Q(t) . Since two carriers are 
present, holes and electrons, the total charge Q(t) is composed of two 
parts; hence. 


Q(t) = Q it) + « n (t) 
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Figure 4.9 The upper and lower figures indicate the drift velocity in 
silicon for electrons and holes, respectively. The temper- 
atures in the upper figure are in Kelvin and the lower profile 
is at room temperature. The upper figure is by Taueh, [1958] 
and the lower figure by Dodge, et at., [1964]. 
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Now summing, as indicated above, gives 


X.(t) - X . 

f jj\ C l Ol 

- l q.' 1 

i 1 


where X . is the location where the i*" ionized charge is generated, X. 

CV v 

is the position in which the charge moves away from X Q ^ and 


. = _£ &L 


a 1 oj . 


where q is the charge unit (1.6 x 10 coulomb) and E is the average 
energy required for ionization of a hole-electron pair. Hence in the 
limit , 


-i IM- 


y*>- t M — 

r C K J 


*(t) - x. 


(4.6) 


Finally, J(t) may be obtained from Figure 4.9 by solving the dif- 
ferential equation 


v d -%' % [.(*)] *) 


with initial condition X(o) = L. 


Similar expressions hold for Q n using electron associated mobilities, 
a minus sign for the negative charge and the X(t) initial condition 
= 0 at x = 0. The two charge components have been calculated by 
Quaranta et at. [1965] cjid are reproduced in Figure 4.10 as a function 
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Figure 4.10 The upper and lower figures indicate the charge pulse waveforms i 
a surface barrier detector for electrons and holes, respectively. 
(Data replotted from Quaranta et at. [1965] )„ 
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of time for a type II detector. The calculations were made using a- 

18 -3 

particles with donor concentration (N^) equal to 3.2 x 10 m and 
L = 290 ym at a 200 volt bias. 

The numbering of curves can indicate either various incident ener- 
gies or angles of incidence. Curve numbers 1, 2, 3, 4, and 5 are 
associated with energies 25, 24,7, 23.8, 22.5, and 20.6 MeV, respectively 
for perpendicular radiation or impact angles of 0°, 15°, 35°, and 
45°, respectively for 25 MoV radiation. 

4.4.3 Dead zone region. Two principal factors influence the above 
results; the dead zone layer and the collection efficiency. First, the 
effective dead zone may be determined as follows. By assuming that the 
dead zone is a slab contained within the gold or n * (p + ) surface, 

[Euing, 1962] different travel distances are expected through the dead 
zone depending on the particular angle of incidence [Williams and 
Webb, 1962]. 

If the slab (window) thickness is w (Figure 4.6) then the maximum 
energy deposited [MED) is proportional to the initial particle energy 
minus that which is lost in the window; hence 

dE 

mD s E 0 ~ W ~Zx 


By varying the angle of incidence 0, it becomes 


MED 


dE 

w o 

COS0 ~3x 


(4.7) 


where 0 is measured from the normal (Figure 4.6). 
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The differential energy for angles compared to the perpendicular 
then becomes 


A E - W 


dE 

o 

dx 


1 

cose 


(4.8) 


Therefore, by illuminating a collimated monoenergetic beam of 

dE 

particles upon the surface, the value of w can be found by measuring 
the angular dependence (a slope) of the peak pulse height since it is 
proportional to MED. 

Ray and Barnett [1969] utilizing this technique have determined 

an energy loss of 2.59 keV for a 10 keV preton beam incident upon a 
_2 

40ygm cm gold surface barrier detector. 

4.4.4 Collection efficiency. The second loss mechanism is attri- 
buted to the fact that not all of the liberated charge from an ioniza- 
tion source is capable of continued motion in the electric field across 
the depletion region for some recombines. The collection efficiency, n, 
is defined to be equal to Q , the collected charge, divided by Q r , th ‘ 
charge liberated within the depletion region. An upper limit to the 
charge collection time may be established if the depletion width L and 
the appropriate carrier velocity (Figure 4.9) are known. The recombina- 
tion is essentially of three mechanisms: 1) recombination via excita- 

tions; 2) direct radiative recombination [Chynoweth and McKay, 1956]; 
and 3) recombination centers "trapping" [Shockley and Read, 1952], 

An analytical form for collection efficiency for particles with 
range distances short compared t.c depletion depth L is given by 
Miller et al. [1962] as 
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n 



-l 


2y «x *6 
%> t m 


(4.9) 


where y. and x. are the hole or electron mobility and lifetime, respec- 
tively, and is the maximum electric field. From this, one deduces 
that, first, the distance, "trapping length", eyx, should be long com- 
pared to d and, second, that "trapping" losses are independent of 
applied bias. 

Since "trapping" is a nonuniform process, a distortion effect is 
introduced yielding poor resolution and the so-called "multiple peaking" 
[Miller et al . , 1962] which is the occurrence of two adjacent energy 
peaks observed from a monoenergetic beam. Additionally, multiple peak- 
ing can be produced by charge multiplication caused by carrier injection 
at the device contacts. Figure 4.11 illustrates this multiple peaking 
behavior for the detector displayed in Figure 4.1. Both photographs 
are taken from the multichannel pulse height analyzer for 60 keV ener- 
getic electrons with the electronics described in Chapter 5. The lower 
photograph is the logarithm of the top linear plot. The peak to the 
far left is the noise. 

In summary, the depletion region electric field is, on the average, 
linear in type I and II detectors (Figure 4.6); however, because of fab- 
rication inhomogeneities, perturbations exist causing the ionized carriers 
to traverse a somewhat larger path and, consequently, have less chance of 
survival for a fixed lifetime, x. The fractional loss of carriers by 
trapping is given by Shockley and Read [1952] to be 



n t 


where x is the recombination lifetime and t 

P 


is the plasma time. 


(4.10) 
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Figure 4.11 The u{ i ind lower 1 are a linear md log 

plot, respectively, of the sol id -state detector 
output analyzed on a 400-channel pulse height 
analyzer. The multiple peaking effect is observed 
at the right, the high peak represents the 60 keV 
incident electrons and the far left peak is the 
noise buildup. 
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4.4.5 Plasma time and vise time. The plasma time is defined to be 
the period that the particle track is held together by its own electro- 
static forces. If the electric field is sufficiently loi. it is Table to 
fully penetrate the high density plasma track. Hence, the charge collec- 
tion time is lengthened. 

Tore and Seibt [1967] have resolved the following functional rela- 
tionship for t in a coaxial detector. 


where 


an r. r 

t = - ---- - -U- 

v 2 v 

‘ yee 

0 


(4.11) 


- 19 

q is the electronic charge (1.6 x 10 coulomb) 

r. is the relative dielectric constant (12 for silicon) 

- 12 - 1 

e is the vacuum dielectric constant (e = 8.85 10 Fm ) 

o o 

is the radius of the original track 
r is the radius over which the track charge is distributed 
by unit transit time 

Meyer and Langmar [1965] have resolved the following experimental 
relationship for a-particles in a surface barrier detector. 


a 

2.2 
z o t 


where a * 1.9 x 10 ^ + 15% sec cm ^ and e = 
The time required for the collected charge 


its maximum value is called the rise time, t . 

’ r 


(4.12) 
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Figure 4.12 Illustrated is the charge pu*. isetime T r as a function of either R/X 0 ratio or 

impact angle, 6, for various impurity concentrations U ^ in a surface barrier detector. 
(Replotted from Quaranta et al . [1965]). 
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[Quaranta et at. , 1965] compares rise time as a function of impact angle 

ft 

6 or ^ ratio for differing impurity concentrations N^, in a surface 

barrier detector. An intuitive understanding for why the rise time falls 

E 

rapidly for increasing angles or small j- ratios is that the ionization is 
deposited closer to the boundaries of the depletion region. Since there 
is less distance to traverse ( R sin6 for nonperpendicular radiation) the 
collection time would be less giving a faster rise time. More funda- 
mentally, the rise time is composed of three components in a time field 
yielding 


U R 


. 2 .2 * 2 
t v + *p + t d 


1/2 


(4.13) 


where t = rise time of the charge pulse due to transit time only. 

t = plasma time. 

F 

tg = rise time associated with capacitance and resistance of the 
detector equivalent circuit. 

4.4.6 Channeling. Crystallographic symmetry with respect to the 

incoming radiation, causes differing quantum efficiencies of ionization 

dE 

resulting in a nonuniform deposition rate. The most symmetrical cut 
of a crystal is normal to the axis of the crystal growth, the [111] direc- 
tion, thus giving the highest percentage of channeling. Gibson [1966] 
has transmitted 4.9 MeV protons through 35 ymSi crystal at directions 
parallel to the [111], [110], and "no symmetry" directions (Figure 4.13). 
Deamaley [1964] has determined that wafers should be cut at 10 to 20 
degrees from the [111] plane in order to minimize symmetry (channeling). 
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Figure 4.13 Channeling io an important factor in the orientation of a 
detector crystal and particle calibrations for wide angle 
applications. The above curves show the strong dependence 
on the relative angle of incident radiation. Data taken 
from Gibson [1966] . 
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4.4.7 Detector equivalent circuit. The output charge generated in 
the depletion region is modified by the surrounding conductors, contact 
points and internal electrical semiconductor properties, i.e. polariza- 
tion. Simplifying, the equivalent circuit may be thought of as containing 
four impedance components associated with the four distinct regions of 
the detector; the depletion region, undepleted semiconductor region, con- 
tact points, and the external connectors and shielding. A more rigorous 
treatment may be made employing transmission line equations, especially 
for the high frequency case. Each section, because it is capable of 
charge storage across a potential, constitutes a capacitance along with 
the ohmic loss. Figure 4.14 illustrates the equivalent circuit where 
C q and C g are evaluated in the same way as that of a plane parallel 
capacitor. 


0^ = — — the depletion layer capacitance 


C 
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D-L 

p ~r 

pte Q 


undepleted zone resistance 


undepleted zone capacitance 


(4.14) 


where p is the semiconductor resistivity. The depletion region resistance 
is negligibly small for a reverse biased diode. The impedance, Z (contact 
point), can considerably influence the pulse and shows considerable 
distortion at lower temperatures. The connectors are primarily capaci- 
tive in nature. Techniques developed to measure circuit parameters are 
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discussed by Quaranta et al . [1966] and Quaranta [1969], A nomogram is 
reproduced in Figure 4.15 to facilitate understanding of depletion region 
thickness, capacitance, resistivity, and impurity concentrations with 
applied bias. Figure 4.16 illustrates the contribution of electrons 
(a) and holes (b) to the final current pulse for a 400 urn thick surface 
barrier detector using 6000 ft cm n-type silicon [Quaranta et al. , 1967]. 

Figure 4.17 is a typical output current pulse from a surface barrier 
detector for three bias voltages (90, 150, and 250 volts) with amplifi- 
cation utilizing a 3 ns rise time amplifier 

4.4.8 Sensitivity of detector operation to external factors. Ex- 
ternal effects of interest include 1) radiation damage, 2) light sensi- 
tivity, 3) temperature, 4) humidity and contaminants, 5) pressure changes, 
6) mechanical shock, and 7) magnetic fields. 

Radiation, if extreme, is harmful to detectors since nuclear modi- 
fication results and crystal structure is perturbed by movement of atoms 
into intersitital positions. In rocket-borne detectors ionospheric 
radiation conditions are moderate and short; operation psriods result 
in negligible degradation. Year long satellite programs through the 
radiation belts record only minor performance loss because of radiation 
damage. 

Light quanta with energies greater than the forbidden band gap of 
silicon are capable of exciting charge carriers ana thereby increase the 
reverse current of the diode and ultimately the noise. Figure 4.18 shows 
the noise introduced by a light source on the detector near the breakdown 
region. Surface barrier detectors are the most sensitive to light 
because of the thin gold layer; the thicker the layer the less light 
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Figure 4.16 Current pulses due to electrons (upper) and holes (lower) in a 400 ym thick surface barrier 
totally depleted detector. Horizontal sensitivity is 1.75 ns/main division lower. Vertical 
sensitivity is the same. Data replotted from Quaranta et at . [1967]. 











Figure 4,18 Preamplifier output illustrating the breakdown phenomenon 
resulting from a bright light source (eg. the sun). The 
vertical scale is 0.1 volts/division and the horizontal 
is 2 msec/division. 
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sensitive. Special coatings can be placed over the gold to reduce light 
passage. Greater bias voltages can also be applied to a detector in the 
dark. Figure 4.19 is a plot of noise verses light intensity for the 
detector shown in Figure 4.1. 

Temperature effects will be discussed in Chapter 5 in greater detail. 
The noise is gaussian and varies as the expotential of the inverse of 
temperature . 

Humidity and vapors of organic solvents are able to react with the 
silicon (e.g. silicon hydride) or gold surface producing harmful effects. 
Silicon detectors should be stored with some drying agent in a sealed air- 
tight container. Whc*. detectors are pumped down in vacuum systems, mer- 
cury diffusion pumps should be avoided since mercury reacts with gold. 

Surface leakage currents can be a dominant noise source as compared 
to space- charge generated currents in the presence of water vapor or other 
contaminants. Figure 4.20 shows the marked effects of dry oxygen and wet 
nitrogen on the surface leakage current of a diffused n + p junction 
[Buck, 1961]. Accumulation of moisture and/or ionizable impurities near 
the surface can break down the intense electric field region in micro- 
plasma. If a detector is stored over a period of time or is exposed to a 
contaminant the bias voltage should be raised slowly allowing the con- 
taminants to drift out of the field region. Figure 4.21 illustrates the 
intermittent behavior of microplasma breakdown. The bottom trace is the 
preamp output showing saturation of the signal over most of the trace. 

The top trace is the postamplifier output again showing similar 


behavior. 



Channel I Counts (Ramps- sec 








Figure 4.20 Reverse current versus reversed bias voltage showing the 

strong dependence upon the external air environment. Data 
taken by Buck [1961]. 
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Figure 4.21 Microplasma breakdown is shown in the alternate trace of 
the preamp output (.top trace is at two volts/box) and 
postamplifier output (bottom trace at 5 volts/box). 

During short time intervals the detector is seen to recover 
from saturation noise to its normal operation. 
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Mechanically the detectors are quite rugged and show insignificant 
pressure dependence. Their mechanical stability is similar to 
transistors. 

The solid-state detector is independent of magnetic fields of sev- 
eral thousand gauss because of the intense electric fi?ld strength dis- 
tributed over the very small depletion region. 
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5. DESIGN OF A SOLID-STATE PARTICLE DETECTOR 

The approach will be taken, as in the actual design, of first setting 
down some of the essential rules and design philosophy for an overall 
low noise system moving then into the working block diagram of the detec- 
tion process. A simplifying FORTRAN computer program (NOISE) may be 
employed at this point to aid in the design for a particular detector, 
preamp, filter, and postamp operation. After establishing the desired 
parameters, a description will be given the chosen design. Finally, 
some attention will be given to future improvement. 

5.1 Noise Generation and Reduction in Semiconductor Detectors 

An outstanding problem in low energy particle detection is the 
reduction of noise. Not only does minimizing noise enhance the low 
energy threshold of the detector but it increases the overall resolving 
power for all energies. 

5.1.1 Thermal, shot, ccid leakage noise . Referring to the detector 
equivalent circuit of Figure 4,14. there are essentially three generators 

of noise: the thermal noise of the iesistors R and R , shot noise 

p s 

(also called space charge noise) and noise associated with the leakage 
current I The resistance, R 9 is usually very large and its thermal 
noise contribution can be neglected for the following reason. 

From basic statistical arguments, the mean squared noise voltage is 
given by the relation 


v 2 = m r Lf 


(?.i) 


i 
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where k is the Boltzmann constant, T is the absolute temperature, R the 
resistance, and A/* the frequency range over which the measurement is 
made. Now, applying Kirchoff's voltage law to a resistor (with equiva- 
lent series noise generator) shunted across a capacitance, C, the value 
~ 2 

of y 4 ,,„ across the parallel combination is given by 


l _ 4 K T R AF 
NR 1 + (tuflC) 2 


(5.2) 


and the total mean square thermal noise is obtained by integrating equa- 
tion (5.2) over all frequency intervals. 


4 K T R d f _ KT_ 
1 + (2tt f R C) 2 C 


(5.3) 


Equation (5.2) clearly shows that as R increases V ^ decreases 
providing C is not too small. 

The series resistance is likewise shunted across a capacitance 
2 

and the noise, v is usually negligible even though R^ is much smaller 

2 ~ 2 
than Rp . The reason why the V ^ is comparable or larger than the V ^ 

for the parallel branch is that the series capacitance is much smaller. 

Several components contribute to the leakage current. First, 

leakage around the periphery of the detector in the high electric field 

region produces an increased level of low frequency noise. Imposing 

a guard-ring and suitable epoxy around the edges will considerably reduce 

this leakage. No adequate model has been developed to explain this source. 

The second major contributor to the leakage current is shot noise 

flerived from random fluctuations in the number of charged carriers emitted 
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.from a surface. By applying statistical mechanics again, the functional 
dependence is obtained 

I s 2 = 21 L q A f (5.4) 

where I is the detector leakage current and A/ is the frequency interval 

in which the measurements are made. 

By similar analysis as made for equation (5.2) employing Kirchoffs 
2 

voltage law the v ^ across the parallel combination of the noisy current 
source and equivalent capacitance is 
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(5.5) 


and the total mean square thermal noise is obtained by integrating 
equation (5.5) over the appropriate frequency intervals. 
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Figure 5.1 gives the noisy equivalent circuit. Figure 4.19 shows 
the leakage contributions of sho^ noise and peripheral leakage due to 
contaminants . 

5.1.2 Detector bias regulation . Over the life of the detector the 
leakage current is expected to vary and thus precautionary measures 
should be made to regulate the applied bias voltage. Figure 5.2 is a 
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Figure 5.2 Circuit diagram employed to regulate the bias current for a surfa~e 
barrier detector. Factors known to significantly change the current 
are temperature, radiation, and long-range storage. Proper bias 
insures correct preamp matching and minimizes noise (Circiu.: taken 
from Bakke et al. [1969]). 
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detailed circuit diagram employed by Bakke et al. [1969] to regulate a 
surface barrier detector and would be an excellent improvement over the 
present system, to be described later. 

The depletion region width (and therefore capacitance) is directly 
related to the applied bias (equation 4.16) and consequently to the 
thermal noise. Figure 5.3 is a plot of noise versus applied uxas oicage 
for the Ortec detector displayed in Figure 4.1 

To measure the bias voltage accurately, an extremely high input 
impedance meter must be employed since the detector resistance is 

g 

approximately 5 x 10 ohms. The circuit (Figure 5.4) uses an operational 
amplifier with an input resistor of 2 x 10^ ohms. 

5.1.3 Mathematical treatment of noise and its temperature varia- 
tion. The most significant factor in reducing noise is to lower the tem- 
perature. To facilitate a fundamental understanding „f the arguments 
involved, an investigation will be made of the spectrum characteristics 
of noise. 

Since noise in a material originates from the unequal balance of 
random mutually-exclusive charge-carrier events, the guassian distribu- 
tion applies. It is expressed in mathematical form as 


P[E)dE - — exp 

ITO 



dE 


(5.7) 


where P(E) is the probability that in an energy interval E to E + dE 
N. 

there exists — charge movements across the band gap E , AL is the charge 
o th 

movements in the i interval and N is the maximum charge movements made 
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Bias Voltage (Volts) 


Figure 5.3 Surface barrier detector total noise versus bias 
voltage. Designed bias is 130 volts. 
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in the interval, a is the standard deviation, and y is the offset of the 
gaussian peak due to unequal probabilities of generating ( g ) and recom- 
bination (r) of charge carriers. The distribution is depicted in 
Figure 5.5 along with the defining parameters and may be compared to the 
actual distribution taken from the detector utilizing a pulse height 
analyzer. 

The essential statistical expressions for y and o are given by 


= C g-r)N E r 


Ngr E 


(5.8) 


Here E is the mean energy required to generate one electron hole pair, 
y 

For pure noise, i.e. random noise, g = v and therefore y = 0 and 


° - *gJh" E i 


(5.9) 


where q is the elementary charge unit 1.6 x 10~ coulomb and Q ^ is 

equivalent noise generated at the amplifier input. The second equality 

intuitively comes from the fact that the number of carries which jump 

the energy gap E (a capacitance, have a dependence on energy 
2 

equal to Q The next two important relationships for a follows 

directly from their defining properties. 
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60 keV electrons 
FWHM = 9.5 keV 


1 


I 



I 


Figure 5.5 Experimentally-derived distribution for 60 keV electrons obtained from a 
pulse height analyzer. Each channel represents an interval of 0.317 keV 
giving a FWHM of 95 keV . In this case p is channel 198 representing 60 keV 
and a = FWMH/2.35 = 4.05 keV. 
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Solving for o one obtains approximately 19 mV which agrees close 
with 19,5 mV measured. The small error is due to incomplete specification 
of the complex input impedance. 

Now the first channel discriminator is set near the tail of the 
gaussian distribution and therefore o is of considers le importance, 
especially since a = /( T) . To start with, the probability of an event 
occurring greater than some threshold energy, E^, of the detector must be 
determined. The result is obtained by integrating over P{E)dE distribu- 
tion from E t to “ or 
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The integral values for the $(a) are tabulated and found in any 
mathematical handbook. For the designed detector, the threshold energy 
was set at 11.6 keV hence z - 11.6/3.6 * 3.23 yielding fTom the tables 
$(3.23) = 0.99936 or P(> E f ) = 6.4 x 10' 4 . Therefore, the total number 
of counts from equation (5.9) is 9.4 x 10^ requiring the number of 
counts greater than E ^ to be N ^ = 6000 counts per second. This agrees 
with the experimental results displayed in Figure 5.6 for a temperature 
of 30°C. 

/ 2 

To deduce theoretically the / V vs. channel one count rate, 
equations (5.9), (5.10), and (5.14) must be used. For a change in the 
root mean squared voltage, there is a corresponding linear charge in 
the standard deviation (equation 5.10). The number of counts for 
energies greater than E , W( E ^) , is thus given by 


H> E t ) = P(> e t )n t 




(5.15) 


Substituting for a from equation (5.9) and (5.12) gives 

E, 


n> e t ) = 


1 - erf 
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(5.16) 


Now using i? = 14 mV, equation (5.1b) gives a count rate of 60 
counts per second which agrees beautifully with the experimental curve. 

The temperature dependent curve of Figure 5.6 is likewise related 
to the number of counts above E~ by equation (5.15); however, a is now a 
function of temperature. The particular functional form will be 
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NOISE mV rms 



Figure 5.6 Total noise counts recorded above the lowest channel threshold 
voltage, 0.83 volts for various temperatures and root mean 
square voltages. 
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investigated after a study is made of the preamplifier noise contribu- 
tion. From Figure 5.6 it is apparent that o(T ) must be some function 
which is less sensitive to T than the linear dependence of equation 15.10). 
Referrin , to equations (5.18), (5.9), and (5.15) the general temperature 
dependence of the channel counts above a threshold energy can be 
formulated. 

5.2 Preamplifier Noise and Detector Poetamplifier Interface 

3.2.1 System amplification philosophy. A complete block diagram of 
the detection and amplifying circuits is represented in Figure 5.7. The 
current introduced into the preamp from an energetic particle is inte- 
grated in the first stage by applying capacitive feedback. The result- 
ing charge pulse is shaped and amplified in succeeding stages. The 
preamplifier is a double differentiation and double integration pulse 
shaping network with a low noise high input impedance field effect 
transistor employed in the first stage, Postamp] ification is made for 
final pulse height discrimination. The Norton equivalent circuit is 
substituted for the detector. 

Experiments show that the Norton admittance is primarily capacitive 
in nature where upon simple analysis yields the first stage output. By 
making a Thevenin equivalent to the left of the first stage and assuming 
tne coupling capacitance and resistive feedback are chosen to have minor 
influence on the output, the basic relation for operational amplifiers 
applies; namely 
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where and B ^ are the Thevenin equivalents, given by 


if + 

e t = ~ — : — 

‘ n- + 3 WC 

r b 


R, 


Rr , JWC 
_ B _ egv 

T 1 

R B + juC 


This can be reduced to the following relation 


E 

o 



I[t)dt = 



(5. 18) 


The collected charge produces a voltage signal of magnitude Q/ Cy. 

This is independent of the detector capacitance and thus applied bias 

giving improved overall stability. The available signal energy supplied 

2 

from the detector is Q /2 C requiring C ^ to be minimized for increased 
signal-to-noise ratio. 

5.2.2 Equivalent preamplifier noise generator . The field effect 
transistor (FET) used as the first stage of the preamplifier was chosen 
for its inherent low noise performance, simplicity, ruggedness, and 
ability for substantially reducing noise upon cooling. Later improve- 
ments would have the FET mounted adjacent to the detector and both 
cooled, for excmple by dry ice. On an energy basis the sensitivity 
could be improved by a factor of two. The present design utilizes a 
detector designed specifically for rocket-borne applications by 
Nucleometrics and is pictured in Figure 5.8 and in schematic form in 
Figure 5.9. Table 5.1 list the preamplifier specification. 

The proper interface of the detector with the preamplifier requires 
optimization of the preamplifiers sensitivity, time constants, and input 
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„ ure 5 8 preamplifier assembly ami electronics developed by Nuclemte 
specifically designed For detector pictured in Figure l.l. 
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Table 5.1 

Charge Sensitive Amplifier 
Specifications 

Sensitivity: 22 to 110 v/pC 

Time constant: 0.5 to 2 micro-sec 

Load resistor: 100 meg ohm 

Dynamic input capacitance: 5000 pf 

Temperature stability -40°c to +60°c: ±1% 

Gain variation with supply voltages: 1%/v 
Dynamic output impedance: <20 ohm 
Integral linearity to 5.0 v output: ±0.2% 
to 2.5 v output: ±0.05% 

Overload recover x 10 : 30 micro-sec 

x 100 : 50 micro-sec 

x 1000: 100 micro-sec 

For 40 pf input capacitance: Noise: <5 keV FWHM (Silicon) 

Rise time: <20 nano-sec 

With 100 meg ohm load resistor: Noise: <2 keV FWHM (Silicon) 

Power required: Positive: +9 to +12 v, 8 ma 
Negative: -9 to -12 v, 2 ma 
Detector bias: 250 v (max) 

Input/output polarity: Detector input: negative 

Test pulse input: negative 

Amplifier input: positive 

Dimensions: 4.375 x 2.375 x 0.820 in 

Weight: <6 oz. 

Connectors: Input/output (each): Microdot S-93 screw receptacle 

Power /Test/3ias : Cannon DEM-9P 

Environment: Appropriate for use in Nike Apache rocket payloads. 
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impedance. A rigorous analysis requires the aid of computer program 
routines and will be discussed in Section 5.2.4. Briefly, the mean 
equivalent input noise charge, characterizing the preamplifier 
detector is the aggregate noise of the component systems. The two pri- 
mary types of noise sources have been introduced in previous sections 
and are represented by equations (5.1), (5.2), and (5.5). The total 
noise is then the sum of the energy contributions of each of the parts 
being related to the voltages as V 2 ^. The total noise is consequently 
the square root of the sum of the squares of the individual hence 


AkT R At 

Q s C ' + 

N eqv t 2 


4 kT t . 2q ( I L +I c? X 


fl eqv 


(5.18) 


Furthermore, by making the reasonable approximation that the flicker 
noise and the R ^ noise are minor the following simplifications result 
at a temperature of 273°K. 


«,= 1SX 10- 21 *gZ 2 - ^ . 3.2 x 10 - 19 (J +1 ) 
N t K g L 


(5.19) 


The optimum time constant, t for minimizing the noise is obtained 
by differentiating equation (5.19) with respect to x and setting the 
derivative to zero. The equivalent input noise for optimum time constant 
is thus 


x = 0.22C y- ? 9 ? r 

opt eqv I + Ij. 


(5.20) 
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Again it is r>m)arent that decreasing the input capacitance and 
leakage current is desirable along with increasing the trans conductance, 
g, since s 0.7 /g [ Bilger , 1966] for a FET. The optimum input noise 
level is consequently 


N opt 


= 1.3 x 10 


-19 


: r (i + i r ) 

eqv eqv K g 


(5.21) 


The particular equivalent noise input for the preamplifier for 
various input capacitance and time constants 0.5, 1.0, and 2.0 ysec are 
shown in Figure 5.10. The electronics noise contribution can be ex- 
pressed in charge units since this is what the amplifier is designed to 
measure. Either the units are expressed in equivalent electron charge 
(coulombs) or in terms of their spread in a gaussian distribution 
associated with the probability that a monoenergetic particle will pro- 
duce carriers. The latter method uses the full spectral width at half 
the maximum (FWHM) of the distribution. The output noise can alter- 
natively be recorded as a root mean square voltage or current and is 
related to the equivalent charge by equations (5.9) and (5.10). Since 
£ is a function of the material used the FWHM must also be associated 


with a material e.g., FWHM (silicon). 

Two methods employed to experimentally measure the FWHM are by 
direct measurement of the spectral width from a multichannel pulse 
height analyzer or indirectly using the root mean square voltmeter 
method. The latter method is based on applying a generator voltage 

pulse of amplitude 7 through an injection capacitor C, to the 

9 en 




Figure 5.10 Amplifier noise versus input capacitance for preamplifier 

time constants of 0.5, l.C, and 2.0 Msec. From Nucleometrics 
specification sheet. 
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detector preamp input. The signal input charge to the FET is the series 

combination of C. and C 

1 eqv 


^ C eqv + *1 


(5.21) 


and if C, is chosen small with respect to C & the above equation becomes 


Q . s c, V 
8vg 1 gen 


If Q n is the equivalent noise charge at the amplifier input required 
to produce y noise voltage at the output than for an overall calibra- 
tion gain G the following expression is valid 
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and the final result obtained is 


Q . v V C v 
Q si# jms gen ms 

V N ~ V . + V . 

Btg Ptg 


(5.22) 


Referring to the schematic diagram of Figure 5.9, the time constants 
are governed by C 30 , C^, C 32 » and Table 5.2 gives the 

suggested values for selected time constants of 0.5, 1, and 2 msec. 

For the Ortec ruggedized surface barrier partially depleted defec- 
tor (Figure 4.1) the capacitance is 52 pf, the noise is 5 keV FWHM and 
reverse current 0.23 uamps. The of the detector preamplifier is 

_3 

approximately equal to or less than 10 ohms. Substituting these para- 
meters into equation (5.20) yields a design time constant of 0.5 ysec. 
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Table 5.2 

Amplifier time constant adjustment. 


Time Constant (usee) 

n 

CkI 

o 

C 31 

C 32 

C 33 

R 33 

0.5 

100 

3.0 

100 

3.0 

5.1 Meg 

1 

200 

6.0 

200 

6.0 

2.4 Meg 

2 

390 

13 

390 

13 

1.3 Meg 




The gain resistors 


R 32 

R 34 R 35 

R 37 P 58 

5.1 

150 5.1 

0 150 




1 


5.2.3 Pulse shaping network 

The step response to a positive 10 mV test calibration pulse produces 
the output response shown in Figure 5.11 for a sweep rate of 0.5 msec 
per cm and vertical sensitivity of 0.2V cm. The y^^ output noise vis 
measured to be 0.19 volts yielding an equivalent input noise of 
5 x 10 ^ coulombs by application of equation (5.22). For comparison 
purposes, a 10 keV particle loses approximately 3.5 eV per released 
charge carrier (as previously discussed) contributing a total charge of 
4.7 x 10 ^ coulombs ; comparable with the noise as anticipated. 

The approximation made in equation (5.18) that the bandwidth is tw 
is about 1 /t where t is the time constant of the shaping circuits can be 
more fundamentally formulated in terms of the system frequency response. 
Proceeding with equation (5.17) with the stipulation that the flicker 
noise component is of minor significance the translated equality 
becomes 
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Figure 5.11 Preamplifier step response (top trace) for a positive 
10 mV test calibrate input pulse [bottom trace). The 
horizontal sweep rate is 0.5 pscc per main division 
and the vertical sensitivity is 0.2 volts per main 
division. 
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To obtain the overall spectral noise characteristics an integration 
must be perfomed over all frequency intervals Aw of the product of the 
noise contribution times the attenuation index, f(w)/F(o), for the fre- 
quency w where /(w) is the frequency response of the system and F(<?) is 
the normalization factor equal to the maximum value of /(w) or equiva- 
lently the maximum amplitude of the system step response (Figure 5.10). 
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(5.24) 
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Listed in Table 5.3 [Tsukuda, 1964] are the significant pulse 
shaping networks and their associated parameters. The optimum equ. va- 
lent noise charge, (ENC) , is the Q ^ divided by the unit charge 1.6 x 10" 
coulomb. Here it is seen that the delay line clipping (DL) networks per- 
form better than the RC types; however, to construct filters possessing 
variable time constants, as is mandatory in delay line clipping, is not a 
simple undertaking. Improvements here could reduce noise by a» much as 
40% and futile designs should incorporate the delay line clipping methods 
as they become available. 

For the system designed, a RC filter was chosen using double dif- 
ferentiation and double integration stages. Referring to the schematic 
Figure 5.9 the first and second differentiation is governed by R^-,, 

6’ 30 , and R^> R 36 , R^, C^ 2 > and the first and second integration by F J4 , 
C and g» ^ 33 > respectively. The input pole zero cancellation by 

^ 42 * C 32’ 

The frequency response f(w), is obtained from basic circuit analy- 
sis although very tedious in procedure for a complete analysis (NOISE 
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Table 5.3 

Optimum signal to noise ratio obtainable with different pulse shaping networks. Data compiled by Tsukuda [1964J 
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computer solution. Section 5.2.4) the primary response may be made 
referring to the simplified circuit illustrated in Figure 5.7. Trans- 
forming into the Laplace transform space, the output response is qui ce 
simply the system transfer function, H{s), multiplied by the inpu 4 
excitation. 

yS) = H'.a) ya) (5.25) 

The frequency response is obtained from the complex Fourier trans- 
form space which may be obtained by substitution of jw for s in 
equation (5.25). 

The magnitude, |//(ju) | = H[w) , is the corresponding frequency re- 
sponse and tht angle 0 of the complex point tf(ju) is Lhe phase change of 
the output to the input excitatio” 

The step (V. = — ) response of the preamplifier is shown in Figure 

'VS 

5.10 and may be used to graphically evaluate H{'j) - V (jw) | where 
F^(e) is determined from the Laplace transform definition. 

oo 

yjtf) « V o (8) = j u Jit) e~ St dt (5.2b) 

0 

To gain an understanding of the function .1 dependence of the fre- 
quency response an analytical expression may be obtained utilizing the 
simplified circuit of Figure 5.7. Here the transfer function is the ro- 
duct of the transfer function of each of the similar pulse shaping net- 
works. For an operational amplifier with feedback, Z~, and input, 2., 
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impedance the gain and thus the single transfer function is given by 
and the two stage complex transfer function is the product 

S fl Z f2 /Z H Z i2 


B{8) = 


S R lf C n R 2f C i2 

(1 ♦ S T/1 ) (1 + 5 T/2 ) ( 1 ♦ Sx^) (1 + St i2 ) 


(5.27) 


where x^ = 1/(1?^ C^) is the time constant of the first stage feedback, 
etc. For stability reasons the two feedback and input time constants 
are made identical giving a frequency response of 


H(w) = 


2 2 2 

w rJ c. 

£ f: 

,, 2 2.,, 2 2, 

(1+W Tj.)(l+W T • ) 

I 'L' 


(5.28) 


2 2 

This is comparable with the (RC) - (1?C) shaping network suggested 
from Table 5.2. The frenuency response is plotted in Figure 5.12. 
Substituting equations (5.18) into (5.24) and evaluating, the equivalent 
noise charge may be acquired. As before, differentiation of Q ^ with 
respect to the time constants and setting the expression to zero will give 
a minimum for the optimal time constants, x^ and x^.. Further 

addition of integration stages improves performance with the theoretical 
limit as n -*• =° of a signal-to-noise ratio of 1.12 as compared to 1.00 
for the cusp response filter. Emphasis is thus made to utilize the post- 
amplification stage as an additional integrator to enhance overall noise 
performance. The step response of the postamplifier is illustrated in 
Figure 5.13 for a 0.1 volt, 1 ysec, negative going pulse. The rise time 
observed is 0.45 ysec as compared to the 0.5 microsecond rise time of 
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Frequency (Hertz) 


Figure 5.12 Typical frequency response for a double integration 
double differentiation pulse shaping network re- 
presentive of equation (5.28). 
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the preamplifier. Future optimization is possible by employing capacitive 
feedback in an integration mode. 

The final item of interest is the functional time dependence of 
the output signal. Tiis is resolved immediately since the transfer 
function, f (s) , is already available in equation (5.27). For a step input 
the output signal is 


v o {t) = £ _1 <*(«) V s)} 



_ SC* g£ 2 . 

( (1+St^.) (l+Stj) 


(5.29) 


evaluation of this equation using the method of partial fractions leads 
to a time dependence of 


-t/ if. -t/x^ 

v - A e + B e + Cte + Dte 

o 


(5.30) 


where A, B, C, and D are functions of the circuit passive elements. This 
equation expresses the observed output response depicted in Figure 5.9. 

5.2.4 NOISE computer program. To facilitate a comprehensiv eva- 
luation of a sensor-amplifier shaping-circuit electronic system a des- 
cription of a general noise analysis program developed by Motchenbacher 
and Fitcher [1973] is here investigated with particular emphasis on the 
solid-state detector design. Quick, accurate, low cost and low noise 
design is obtained along with a better understanding of the processes 
involved. The program is able to simulate the amplifier network by 
either individual component interconnections or by overall amplifier 
response curves. The six available outputs are listed in Table 5.4. 
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Table 5.4 

Available outputs in computer program NOISE. 

1* Total equivalent input noise over a band. 

2* Input network frequency response. 

3. Input noise at any one frequency. 

4. Input noise versus frequency. 

5. Total noise at the output. 

6. Total system gain response. 
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The source model for the solid-state detector applicable to the 
computer program is represented in the first section of Figure 5.14 and 
the particular element values are listed in Table 5.5. 

The first stage amplifier model (Figure 5.14, Table 5.5) called 
NOSMOD represents the most crucial gain stage. The noise spectrum of 
Epj and I ^ can be approximated by the piecewise linear relationship 
shown in Figure 5.16. 

The final model (Figure 5.15), AMPLGN, is the main amplifier and 
pulse shaping networks. Any arrangement of the curves shown in Figure 
5.16 is valid to represent the particular pulse shaping frequency 
response. For example, N8 is the rate-of-rise in 6 dB/octave steps of 
the response curve. For N8 = 2 the rise is 12 dB/octave, etc. The end 
points for each segment must also be specified. The derived frequency 
response curve is illustrated in Figure 5.12 and provides typical input 
to AMLLGN for analysis permitting one of the parameters to vary in 
equation (5.27) thereby recovering the most optimum configuration. 
Further references should be made to Motchenbacher and Fitcher [1973] 
for a generalized treatment of noise systems. 

5.3 Pulse Voltage Discrimination and Counting Circuits 

The design approach for interfacing the solid-state detector- 
preamplifier output pulse with the telemetry is based on a conventional 
Geiger counter system design in addition to a pulse height discriminator 
circuit. The system, as designed, was intended to be kept as elementary 
as possible fulfilling the experiment requirements. The two primary 
requirements referred to in C lapter 2 are to ascertain the particle flux 
and energy spectrum of energetic electrons so that particle ionization 
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Table 5.5 

Symbol equivalents used in computer program NOISE. 

R = Rl * amplifier input resistance 

r = R2 * real part of sensor leakage resistance 

R r = H3 = load resistance 

Cj - C 1 = amplifier input capacitance 

C - C2 - cell capacitance 

u 

C * Cl - wiring capacitance 

W 

Lp = LI = shunt inductance 

F = F6 = 1// excess-noise corner 

■u 

= J2 » sensor leakage current 
I and = (41c Ttf/R)^^ 2 = thermal noise of F g and F g 
I ^ = ( 2<7 A - shot noise of cell current 

E and I - E and I = amplifier equivalent noise mechanisms 

n n 
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NOSMOD 



Log frequency 


AMPLGN 



Figure 5.15 Frequency response models for use in subprograms NOSMOD 
and AMPLGN. 
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Figure 5.16 System representation of the solid-state particle detector experiment. The particle 
produces a charge within the detector, whereupon it is amplified and pulse height 
discriminated. Fach pulse within the discriminator range triggers a digital 4 bit 
counter and is then converted into analog form for telemetry transfer. 
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rates coulo be established. To most simply accomplish these investiga- 
tions a six channel energy analyzer circuit was decided upon since pre- 
vious indications [Potemra, 1973] suggested that the spectrum followed a 
monotonically decreasing count rate verses increasing energy. More 
elaborate procedures for data manipulation and reduction are available 
options to the system as familiarity, feasibility, and scientific 
objectives dictate. 

Further concern of the design includes; 1) shock and vibration 
protection, 2) shielding from interferring payload equipment and high 
power radar, 3) ease of replacing and checking units, 4) operation under 
low pressures, 5) temperature insensitivity, and 6) noise immunity. 

A description will be given of the method and hardware implemented 
to apply the above design approach. To properly aid in understanding 
this material an effort is made to include an ample number of photos; 
illustration, and oscilliscope displays. 

First, the electronics will be examined as to system operation, 
noise immunity, design, construction, and future improvements. 

Secondly, the mechanical construction will be discussed in Section 
5.4 and reference should be made here as to component location, 

5.3.1 Counting circuits operation. The fundamental operation of 
the solid-state energetic particle experiment is system represented in 
Figure S.16. Incident particles impinging on the semiconductor detector 
create electron hole pairs proportional to the particle energy. A 
specially designed charge preamp described in Section 5.2 senses this 
current pulse and amplifies it accordingly. The signal is height dis- 
criminated, restoring the particle energy information for subsequent 
particle count and data transmission. 
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| The pulse height circuit uses six series-connected voltage comparators 

as a window discriminator. Each signal is then shaped and accumulated in 
a self-resetting 16-count 4-bit binary counter. This information is then 
converted to analog form for interface with telemetry VCO. 

5.3.2 Noise immunity. One of the most challenging problems 
encountered in the design is elimination of unwanted spurious pickup. 

The three primary noise sources are interferences from the high voltage 
power supply (30 V to 300 V converter) , ground loop modulation, and power 
lead pickup. The final consequence of noise exhibited in the circuit is 
extraneous pulse counts, thereby reducing overall sensitivity and 
accuracy. 

Power supply noise can most simply be reduced by using decoupling 
capacitors. A 0.1 yf disk ceramic capacitor is mandatory for RF fre- 
quency bypassing and a large electrolytic capacitor for transient 
suppression. Further, noise suppression is obtained by bypassing each IC 
with a 0.01 yf disk ceramic capacitor placed as near as possible to the 
IC. Direct radiation pickup from power supply wiring is reduced by 
separating signal and power supply wires and enclosing the power supply 
in a grounded metal shield. 

Ground loops result from the finite resistance of conductors and 
consequently may experience a voltage difference. Possible mechanism 
capable of generating voltage fluctuations are mutual capacitance 
between all bodies, static potentials, changing magnetic flux linkages, 
piezelectric effects, and varying reference plane currents. To minimize 
current loops a ground plane should be used on each circuit board and all 
system grounds tied to a single ground point. For the sensor, preamplifier 
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and electronics the shielding technique should be used as illustrated 
to Figure 5.17. 

5.3.3 Counting circuits design- The detector and associated cir- 
cuits are arranged in the payload on two decks, as shown in Figure 5.18. 
These decks are interfaced through a pin connector to the main payload 
wiring channel, shown in Figure 5.18, for data transfer and low voltage 
(+30V and -30V) battery power sources. 

The upper deck (Figure 5.19) consist of the solid-state detector, 
charge preamplifier, high voltage power supply and voltage dividing 
network. The power converter chops the 30 Volt dc battery pack power at 
a frequency of 5 kHz and steps up the voltage to 300 VDC. The numbers 
in the right-hand vertical box designate the pin connector number. TP1 
and TP2 refer to bias power converter test point and the ca^brate test 
point of the preamplifier, respectively. 

The pcstamplifier design was placed on the lower deck for space 
and power compatibility reasons. Future designs should include the 
postamplifier near to the preamplifier. Also, separation should be made 
of the power converter and regulator on the reverse side of the metal 
deck from the detector and amplifier electronics. 

The lower deck (Figure 5.20) contains three stacked circuit boards. 

The lowest circuit board is used for the power supplies and for the counting 
circuits because the regulators can use the metal deck as a heat sink. 

Three supplies were required; a + 5 Volt for the IC logic., a ± 15 V supply 
for the operational amplifiers and the digitai-t.o-analog converters and 
a + 15 V reference supply for the voltage discriminators. The power 
supply board also contains the postamplifier circuit. The postamplifier 
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ind shielding procedure to minimize ground loops and interfering radiation 
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Figure 5. IS 
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The arrangement of the energetic electron spectrometer in the 
payload of a Nike Apache rocket. The upper deck contains the 
detector, preamplifier and high voltage supply complementing 
the lower deck counting electronics and low voltage power 
supplies. Both decks arc interfaced through pin connectors 
to the main wiring channel shown on the extreme left. 
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' employs a Fairchild yA715 operational amplifier, used in the inverting 

mode. The RC feedback was used for pulse integration in addition to 
pulse gain. Gain response was measured to be linear to 14.5 Volts before 
j saturation was reached. Two compensation capacitors were used to pre- 

vent oscillations, to modify the time constant and to increase stability. 
The basic circuit parameters are discussed in reference materials. 

The upper two circuit boards are identical, each containing three 
complete channels of counting circuits. Further circuit boards may be 
easily added for increasing the total number of pulse height me as ure - 
' ments. The input pulse to the pulse height analyzer circuits is essen- 

tially the inverse of the preamplifier pulse displayed in Figure 5.11. 

I 

I Some modification within the negative portion of the postamplifier output 

occurs (Figure 5.21) due to coupling from the ADC, and shaping circuits. 
This coupling is associated with instabilities in the negative portion 
of the transfer function. However, it is not of concern since the pulse 
height is not affected. No amendable solution was found which did not 
I modify the original pulse. 

.• For voltage discrimination a st<. ' ^d array of voltage comparators 

were employed. Each comparator was re. renced at a different voltage 
from a resistive voltage divider network. To save space, a quad com- 
* parator integration circuit was used on each counting circuit board. 

| Hysteresis was added to reduce the sensitivity of the transition region 

I 

to approximately 3 mV. The typical output from a comparator is illus- 

l 

trated in the top trace of Figure 5.22 for the input pulse in the bottom 
trace. The comparator output is observed to be delayed approximately 
' 12 usee from input pulse trigger. 
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Figure 5.21 Output signal (upper trace) from the postamplifier for 
a 10 inV test calibrate signal (lower trace) supplied to 
the preamplifier. The horizontal sweep rate is at one 
microsecond per main division. The oscillations apparent 
in the negative portion of the pulse are the result of 
coupling from the ADC and shaping circuits. 
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The next stage is used for triggering and predetermining the dead- 

times specified for each channel. When a negative-going signal is 

received from the comparator the one-shot locks on for a set time and is 

not affected by subsequent pulses during the dead-time period. This 

eliminates effects from multiple triggering as well as defining the pulse 

shape, and drives the succeeding counter stage. The value of dead-cime 

chosen is determined from the data capability of the telemetry bandwidth 

when saturation begins. For IRIG telemetry channel 20, approximately 

2 x 10^ counts per second can be transferred accurately. Since the counts 

are grouped into ramps of 16 counts, the maximum resolving power is 
3 4 

16 x 2 x 10 or approximately 3.2 x 10 counts per second. This dictates 

4 -1 

a dead-time of (3.2 x 10 ) seconds = 30 usee. 

4 

For counting rates approaching 3.2 x 10 statistical analysis pre- 
dicts the correction factor for particlc-s which are within the built-in 
dead-time. Figure 5.23 illustrates this behavior. The formulas may be 
referenced in statistical handbooks. The correction factor for a 
dead-time, t^, actual count rate of N and observed counting rate of 



Because of the nature of the particle spectrum, higher energy 
channels are not expected to saturate; therefore, reducing the need for 
longer dead-times for the smaller bandwidth channels. The dead-times are 
also staggered to reduce simultaneous surges in the ADC and counters 





counts registered irui . I.:**, electronics for a random count rate 
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causing distortion effects eluded to in Figure 5.21. Figure 5.24 illus- 
trates the dead-times associated with the four high count channels. The 
sweep rate is 5 ysec per main division where the tcp trace is the highest 
count channel. 

The shaped pulses are next counted in a 4-bit binary counter and 
paralleled across to the ADC. For high counting rates an improvement 
can be made to enhance the counting capability 16 times by employing an 
8-bit (or at least greater than 4-bit) binary counter and resolving only 
the most significant four bits to the ADC. Each ramp is therefore repre- 
sentative of 256 counts giving an overall count rate for channel 20 of 
2 K X 256 = 5.12 x 10 5 sec' 1 . 

The binary count received from the counter is transformed info an 
equivalent analog representation by use of a MC1508L-8 ADC and a quad 
operational amplifier. The operational amplifier is necessary to con- 
vert the current signal from the ADC to a voltage signal of sufficient 
magnitude for interface with telemetry voltage controlled oscillators 
(VCO) . Typical output from the electronics for a random input flux is 
displayed in Figure 5,25. The circuit board layout is pictured in 
Figure 5.26 for the power supply board and counting electronics. 

5,3.4 Future improvements of counting circuits . The present design, 
although satisfactory, has two limitations: 

1) Each spectrum interval requires one IRIG telemetry 
channel (consequently the total number of channels 
available to the detector experiment is limited) . 

Additionally inefficient use is made of telemetry 
space. Fine structure may also play an important 
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Figure 5.24 The above four pulses show the staggering utilized in 

the pulse shaping circuits. The upper pulse corresponds 
to the high count rate channel and is set for the maxi- 
mum bandwidth for 1RIG channel 2U, Horizontal sweep 
rate is 5 usee per main division. 



Figure 5.25 Output signal supplied to telemetry for the high count 
rate channel. Each ramp represents lb counts. The 
counts shown here are the result of noise pulses greater 
than the channel threshold. The horizontal sweep rate is 
10 m sec/ Box yielding a count rate of approximately 
55 counts per second. The temperature is 10°C and the 
rms noise is 14.2 mV. 




Figure 5,26 Meet runic ' j t * e Liiuitinj circuit board (upper picture) 
and power s upp 1 y hoard (.lower picture). 
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role in determination of ionization rates requiring 
a detailed spectrum and ample channels. 

2) Pitch angle data reduction is a quite complicated 
procedure with the present data format. 

These problems are eliminated in an advanced design outlined in 
Figure 5.27 which is capable of 64 spectrum channel coverage. All of 
this information can be integrated on one IRIG channel of an FM tele- 
metry system. The provision is additionally made to have the counted 
data synchronized with the aspect magnetometer for simple pitch angle 
analysis . 

The principle of operation is similar to the system described in 
Section 5.3.1 with the following important changes and additions. The 
input pulse supplied from the postamplifier is tracked and made to 
lock on to the peak of the pulse. An analog- to-digital converter reads 
the steady-stare peak amplitude and resets the peak detector. The digi- 
tal equivalent of the measured peak can then be transferred to a PCM 
telemetry network, if available. However, telemetry space can be signi- 
ficantly reduced if spectrum is internally analyzed. To accomplish this, 
:he 6- bit binary signal is converted into one of 64 bins corresponding to 
the magnitude of the binary word. Each pulse transferred to a bin is 
individually counted. The first and second most significant group of 
4 bits is then converted into two analog voltages, respectively. The two 
analog signals comprise a two digit hexadecimal number resulting in 
128 an;* Log output signals. These 123 signals are multiplexed at a rate 
of 2000 samples per second, the bandwidth for IRIG channel 20. Each 
counter is reset sequentially in step to give a counting period and 
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Figure S,2 7 Auvamed design for telemetry transfer of particle information. Each pulse from the 
postamplifier is height-analyzed and counted in such a manner for compact utilization 
of an IRJG channel. 
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scan period of 128 divided by 2000 = 63 msec. The rocket spin rate is 
approximately 5 rps or 200 msec giving approximately three counter scans 
per rocket rotation. To improve the position sailing the multiplexing 
rate must be increased (i.e. use two IRIG channels), the spectrum infor- 
mation reduced or, to simply reduce the counting period. If the period 
is halved the position information is doubled, etc. The divider network 
is easily adopted to this operation. The basic reasoning in the latter 
approach is that not every particle penetrating the detector need be 
counted but only periodic samples during each revolution. The final 
output would therefore have three position maasurements per second over 
a reduced scanning interval. For spectrum scans of 32 energy channels, 
there would hirty positions measurements per second for pitch angle 
information and 32 spectrum scans per second all of which could be 
transferred as IRIG channel 20. 

The pulse-height analyzer circuit is a non-inverting peak detector 
based on a voltage-following principle. Referring to ie circu. iia- 
gram (Figure 5.20) diodes D1 and D2 are used for transient protection 
and for compensating nonlinearities introduced by D3. Capacitors Cl and 
C2 are used to compensate the loop for fast tracking applications in 
which overshoot is present. The response speed is also dependent on 
C^,the holding or storage capacitor. Since speed is important and 
storage is only needed for a few microseconds, C^ is correspondingly low. 
Reset of C^ to zero voltage is initiated by the FET saturating after a 
predetermined time interval governed by one shot 1. The second one shot 
is set slightly less than the ^irst one shot so that the ADC can read 
the stored voltage. 
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The binary to decimal converter is fabricated as a complete unit 
integrated circuit. The counters and AC are operated in the same manner 
as discussed in Section 5.3.4. The multiplexer is again a standard off 
the shelf item and must be capable of sequentially reading twice the 
number of spectrum channel. The clock is a crystal controlled oscilla- 
tor with an attached counter for dividing down the clock rate so the 
accumulators can be reset. 

5.4 Mechanical Layout of Detector System 

The mechanical design required ir+egrarion of the varrous compo- 
nents in such a way that (1) electrical interference was minimized, 

(2) flexibility of replacing units and easy subsystems checkout could 
be made, (3) that precautionary measures exist to resist shock and vi- 
bration and (4) that temperature and pressure charges are insignificant. 

For a versatile design, the conclusion is reached that the system 
should be constructed on two decks; the upper deck to contain the ampli- 
fier and detector stages and the lower deck to contain the counting and 
height discrimination electronics. The option for interchanging elec- 
tronics packages could then be made as design improved or scientific 
objectives dictated. Both packages communicate by use of pin connectors 
to the main rocket wire tray as previously discussed. The upper and 
lower decks are pictured in Figures 5.28 and 5.29 

To guard against vibration, shock, and pressure each compartment 
is encapsulated. Th^ foam used, Eccofoam FP, is airtight to maintain 
operation during sudden pressure changes. Furthermore, the foam in 
its solidified state rigidly bonds all the individual components to 
prevent vibrational fatigue and absorb shock impulses. 
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Figure 5,28 Upper deck of the energetic spectrometer* The bias power 
supply is located in the can behind the detector mount* 

The preamplifier and pul-e shaping elec onlcs are located 
beneath the platform* 




* 





u counting electronics. The lower board contains 
ircuits . 
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To protect the payload and, most importantly, the sensitive surface 
of the detector during initial flight two doors are situated on the rocket 
shell. At an altitude of approximately 40 km the doors are ejected, 
exposing the detector. Figure 5.30 illustrates the location of the 
detector relative to the door opening. The detector mount is shown in 
Figure 5.31. 

To further increase the resolving power and directional capabili- 
ties of the detector future improvements should include provision for 
cooling the detector and adoption for collimator tubes. The low energy 
threshold and resolving power could be halved, to about 6 keV, if the 
detector were to be cooled by dry ice. The basic problem here is not 
the difficulty of cooling the detector but more so the application of a 
seal to prevent condensation, a harmful contaminant, to the detector 
surface. A possible solution for sealing is to mount the detector in 
an evacuated cylinder, hiving the same diameter as tne detector. 
Immediately after door ejection the front comer of the detector could be 
removed by solenoid action. The chamber is thermally insulated from 
the mounting platform while thermal contact is made with a dry ice 
reservoir surrounding the perifery of the cylinder. 

Various collimator systems have been devised. Basically, aper- 
tures should be positioned as shown in Figure 4.3 to reduce scattering 
off the walls. Brass has been used to construct these apertures 
because of the greater likelihood of particle absorption. Further adap- 
tions include magnets of approximately 1000 gauss to discriminate between 
protons and electron particle content. 
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where C is the total equivalent capacitance referred to the input and 

eqv — 

v vms is the ec l u: *- valent n oise, \Jv generator referred to the input. 

FWF.M is the full width at half the maximum of the distribution. The final 
relationship for o will be derived later; however, it is advantageous to 
state it now as 


o = E 


C 1 NO V qer 


g a V . 

a ^ stg 


(5.11) 


where \Jv„ n is the rms output noise, C . is the series capacitance through 


which a calibrated signal V is applied (Figure 5.9) to the preamplifier 
input producing an output signal, . 

For the designed amplifier the FWHM was measured to be (Figure 5.5) 
8.4 iceV FWHM (silicon) at a temperature of 30°C. The output rms noise. 


Jv on of the preamplifier at this temperature is 20 mV rms. The cali- 
bration gain, V . /V , was measured to be 34.5 (Figure 5.11) and C. , 

6 * stg ger 1 

is, from the schematic (Figure 5.9), 1 pf. This information leads to 

the conclusions from the above equations that, 1) the standard deviation 

is 3.6, from equation (5.10) or from equation (5.11); 2) the number of 

movement events is 9.4 x 10^ sec"* from equation (5.9); and 3) ^ e q V v m8 

is equal to Q ^ which is 5.2 x 10"*^ coulombs. 

From conclusion three and Figures 5.11 and 4.15 for x = 0.5 sec, 

C is 90 pf yielding v = 5.8 x 10 ^ volts and thus the overall 
eqv r 6 rms 

amplifier voltage gain is 3500. This can be cross-checked by observing 
that the calibrated signal is voltage divided down through C ^ and (7 y . 


hence : 
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6. CALIBRATION OF THE SOLID-STATE DETECTOR 

The calibration procedure involves consideration of the detector 
sensitivity to various kinds of particles with respect to the energy 
and to the particular direction at which the particles enter the semi- 
conductor. The primary theoretical discussions for the above variation 
have been dealt with in Chapter 4. The purpose of this chapter is to 
verify this material with experimental results and obtain the direct 
relationship between the input excitation and the output response with- 
out reference to the system physical makeup. The first section is 
concerned with the directional properties of the detector and how the 
modeled angular dependence calculation agrees with the experimentally 
determined dependence. The second section develops methods for energy 
calibrations using radioactive conversion electron sources with known 
energy emissions. The final section presents calibration results 
obtained from an electron and proton accelerator. 

6.1 Detector Geometrical Factor and Angular Resolution 

One of the most important calibration measurements is to ascertain 
the angular response of the detector incident particles. Experimental 
results indicate that the angular pattern is significantly different from 
the theoretical calculations for electron energies greater than 20 keV. 

The associated reason for this discrepancy is reconciled due to scatter- 
ing of the higher energy particles off of the collimator walls. 

The geometrical factor is a useful quantity which relates the counting 
rate of a detector to the flux, assumed isotropic, to which it is exposed. 
Thus 

G.F . i | (6.1) 
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where N is the detector counting rate and $ is the number of particles 

per square centimeter per second per steradian. Figure 6.1 illustrates 

the detector geometry with the detector mounted in the rocket payload. 

Assuming that all electrons colliding into the wall are absorbed, an 

analytical expression may be derived for the geometrical factor. 

Projection of the circular aperture of the collimator on the plane 

of the detector circular area for a particular angle 0 displaces the 

circles centers by t tan 6 as shown in the lower diagram of Figure 6.1. 

For 0 < there is no vignetting whereby the exposed sensitive area 

2 

is ur cos©, the effective area of a disk viewed from the angle 0. If 
6 > e cr ^ t the exposed area is A cos0 where A is the intersection of the 
two circles illustrated in Figure 6.1 by the section lining. To determine 
A the following relationships are valid 

Sector a + Sector 6 = Triangle ag + 1/2 A 


where 


Sector a 


Sector 8 = 


2 2 
vR a R a 

2 it = 1 

-> 

R“ 6 


2 ? 

Triangle ag = 1/2 R cosa sina + 1/2 r cosg sing 


Solving for A gives 


2 2 
n 2 , sin a . 2 fo sin 6 . 

R (a ^ — ) + r (6 2 — ) 


( 6 . 2 ) 


using the law of cosines yields a and 8 as 




(a - a cos <J>) 'Area = 2a 2 (<£ - sin <p cos <p) 

Figure 6.1 The upper figure illustrates the collimator geometry and 
mounting Hardware for the sol A J-stace detector* Th~ 
lower figure shews the construction for calculation of 
the geometrical factor. 
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The geometrical factor in ster is related to the effective area, 

(6.4) 

since the number of particles counted is the integral of the flux over 

eac»i element of solid angle dSi = 2tt sinB dQ weighted by the effective 

area A . In keeping with equation 6.1 the division by A is required to 

account for the flux, <J>, when the vignetting is zero. For 9 < 0 

tot* ent 



r 9 . . '9 2 222 

cri A cosO dtt = 7T r cos 0 2tt sin9 db = tt r sin 9 . (6.5) 


and for 0 > 0 ^ the appropriate expression is 


0 q i i 

maX A cos0 dn = 2 7i de Sine i? 2 (a - S1 ” 2 - ) + r 2 (S - - ia2B )| cose 

6 . J 0 ^ 1 I 

crit ] (6.6) 


Solving the last integral numerically and substitution for the 
appropriate parameters into equation (6.4) results in the determination 
of the geometrical factor. For the case presented in Figure 6.1 the 
geometrical factor is 1.49 ster 

The angular dependence of the detector is calculated from the normal- 
ized effective area A . Figure 6.2 illustrates thi > angular dependence 
for two collimator geometries. The curve for the beam width of 64 degrees 
represents the geometry of Figure 6.1. Also shown is the experimentally 
measured curve (BW = 82 degrees) using 1L5 heV electrons. The results 
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gure 6.2 Angular patterns showing he relative sensitivity of the electron detector. 

The BW of 64° represents tne geometry of Figure 6.1. The pattern of beam- 
wiuth of 82° was measured experimentally for high energy particles where 
scattering was prominent. 
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clearly illustrate the increased sensitivity due to scattering of high 
\ energy particles off the walls. 

The experimental apparatus used in making the first measurement 
is portrayed in Figure 6.J. A radioactive source (Cobalt 57) coated 
by nylon is used to excite elec rons. The electron source is shielded 
except for an opening having a diameter of 3 mm. The vacuum chamber is 
necessary to keep the electrons from being scattered and absorbed in 
the ambient air. The conversion electron source is directed at the 
solid-state detector and rotated about the detector on a radius of 
13 cm. The detector is connected to counting circuits external to the 
system. The beam width deduced in this manner is 82 degrees. The 
discrepancy between this value and the calculated value of 64 degrees 
is discussed later in this section. 

Subsequently the facilities of the Goddard Space Flight Center 
were utilized. The description of the electron accelerator system 
is discussed in Section 6.3.1. For discussion here it is sufficient 
to note that pulse height analyzer spsetrums were made for various angles, 
0 , of a controllable electron beam. The pulse height spectrum obtained 
for 20 and 40 keV electrons is shown in Figure 6.4. The 20 keV spectrum 
to the left is resultant c’ 6 equal to 0, 20, 30, 40, 50, and 60 degrees. 
The theoretical cutoff angle is 42 degrees suggesting that scattering 
off the walls is the only viable explanation for the 40 keV particle 
counting at angles of 50 and 60 degrees. At 20 keV the theoretical 
cutoff is observed. To better qualify this wall scattering energy de- 
pendence, the experimental results are plotted in Tigure 6.5. Suppression 
of wall scattering effects may be accompli. shed by inserting baffles into 
the collimator section, as illustrated in Figure 4.4. 
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radioactive particle source was made to rotate about the sensitive surface ol the 
l id -state detector and in a vacuum chamber to reduce scattering. Angular patterns 
uld then be made as shown in ligure b 
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Figure 6,4 Pulse height spectrum obtained for 20 md JO keV electrons 

using the instrumentation of the GSIl . The JO kcY spec t rums 
{top photo J are for j articles with angles of 0* 20, 30, 4 0 
and 30 degrees and the 40 he\ spectrums for 0, 20, 30, 40, 30 
and 60 degrees. 
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Particle Energy (keV) 

Figure 6.5 Experimental results were used to obtain this relationship of beamwidth variation with 
energy. These results show that scattering off the walls is significant. 
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6.2 Conversion Electron Sources 

One simple method which can be used to calibrate the output voltage 
signal is to introduce a radioactive source of known energy and measure 
the resultant averaged pulse height or even better yet the analyzed spec- 
trum. The radioactive source employed uses Cobalt-57 in contact with a 
mylar surface. Gamma radiation interacts with the mylar and releases 
energetic electrons of nearly identical energy (minus approximately 5 to 
6 keV, lost in overcoming the work function). For Cobalt-57 the gamma 
energies are 122 and 137 keV producing electrons of 116 and 130 keV , respec- 
tively. The output pulse for this case is portrayed in Figure 6.6. The 
two top peaks can be associated with the above-mentioned decay energies. 

Assuming that the dead zone region absorbs 2 to 3 keV and the detec- 
tor is linear, the calibration curve shown in Figure 6.15 is suggested. 

The method is capable of resolving energies within 9%. 

The error may be halved by employing a pulse height analyzer which 
averages over a time interval to give a most probable peak. Figure 6.7 
shows the double peak produced by Cobalt-57. Again, the calibration 
curve is derivable from this information, assuming linear behavior. 

6.3 Facilities Employed to Obtain Final Spectrum and Energy Information 

6.3.1 Electron and proton accelerator system. Final calibrations 

were made using an electron-proton accelerator at the Goddard Space Flight 
Center in Maryland. The accelerator features a high voltage DC power 
supply and is designed to produce electrons or protons at controlled par- 
ticle intensities with energy variable from thermal to 160 keV . The 
generated beam is directed to a vacuum chamber shown in Figure 6.8 in 
which calibration experiments are performed. The positioning of the 
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Figure 6,6 Output pulse from a Cobalt - 57 radio active source. The two 
pulse heights correspond to the decay energies of 116 and 
130 keV electrons. A quick calibration measurement can 
therefore be made. 
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detector is in line of the beam and located directly above the rotation 
axis of the attached platform. The monitor is positioned adjacent to the 
extreme part and may be rotated in front of the beam for comparison pur- 
poses. All necessary electrical connections are fed through vacuum seals 
to monitors outside the vacuum chamber. 

Figure 6.9 displays the electronics used for calibration. Signals 
from the monitor detector (solid state or channeltron) are preamplified 
at the exterior feedthrough of the vacuum chamber and transferred to the 
main amplifier shaping electronics located on the top of the moveable 
rack. The monitor pulses are then counted and averaged over a five 
second interval period (monitor counter is under top shelf) . The solid- 
state monitor is responsive to electron energies down to 20 keV. For 
energies less than this a linear channeltron shown in Figure 4.3 is used. 
The background counts (ionic feedback) due to the residual gas contribu- 
tes 160 counts per second. 

The preamplified signal from the detector being calibrated is 
similarly fed through the vacuum chamber wall to the counting circuits 
and to a unity gain inverting amplifier located on the top shelf. Signals 
from each channel are connected to a counter and paper-tape punch located 
beneath the monitor counter. The inverted preamp signal was joined to an 
HP noise voltmeter (next to counter), oscilliscope, and pulse height 
analyzer (comprising the foreground electronics rack) . Signals to the 
pulse height analyzer are fed into the peak detector and converted into 
digital form (fourth panel from top) . The digital representation is 
stored in a 4000 channel memory (third panel from top) and interfaced with 
the oscilliscope via the display control panel. Provisions are additionally 
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made for digital printout of channel contents (lowest panel) and polaruid 
photographs of displayed information. The instrument placed on top of 
the PHA is a standardized signal for checking the performance of the 
monitor. 

6.3.2 Channel paseband calibration. The data retrieved from the 
detector counting circuits are dependent upon the efficiency of each 

channel to discriminate and pass voltage pulses. The foundation of 
the problem rest in the distortion of a constant pulse height due to 
superposition of random noise fluctuations. Analytically this can be 
represented by the gaussian distribution (equation 5.7) translated such 
that p represents the energy of the initial incident particle. The 
probability, P(E ) , that the final output pulse will be in an energy 
interval dE is again 

P{E ) dE = — — exp[-(£’-y) 2 /2a 2 ] (6.7) 

2tt a 

Now, 'since the voltage comparators are sensitive to voltages greater 
than a set threshold value, lae passed counts will be all those pulses 
whose probability is greater than E^,. This is the integral of equation 
(6.7) from E ^ to infinity. Representation of this is expressed analytically 

in equation (5.15). and experimentally derived by measuring channel counts 
for incident particle beams every 2 keV (Figure 6.10). 

The major concern is that variation in o resulting from variation 
in temperature, changes the measured efficiency of counting. Figure 6.11 
shows this sensitivity. Accurate determination of the channel efficiency 
curves can be made if the temperature, rms noise voltage, or channel one 
noise count rate is monitored and knowledge is available on the dependence. 



CHANNEL EFFICIENCY CURVES 


1 


! 



Figure 6.10 Experimentally derived channel efficiency count curves which are compai ,!e to the 
derived formulation of equation (5 . 3 5) . These results were made at 63FC. 




Figure 6.11 Theoretical temperature : ation of channel efficiency curves. Equation (5.15) was used 

for the calculation. n ults clearly indicate that temperature is an important 

consideration for prope. o of an effective threshold energy. 
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The actual method employed is to record the noise count rate before 
rocket launch and during descent when the rocket is out of the particle 
region and interpolate between these. For most cases the noise change is 
small and affords only a secondary correction. 

The next important consideration is the effective cutoff threshold 
defined as the energy where half \s maximum number of particles are 
passed per energy interval. Clearly for a continuous energy spectrum 
the effective energy threshold is in the center of the efficiency curve. 
In general, the spectrum is not constant over the threshold region and 
the effective threshold energy is determined by the product of the effi- 
ciency curve, C(E), and particle spectrum P(E) . 


# particles per energy 
interval dE passed by 
the discriminator 


\ = C(E) P(E) dE 


( 6 . 8 ) 


Now assuming the spectrum to be a power law of the form E where x 
is a positive constant and the efficiency curve is linear over the trans- 
ition range, (7(E) = (E-E^) , equation (6.8) becomes equal to 

^E' X (E-Ej) (6.9) 

where o is the noise standard deviation and is approximately the energy 
span between E^ and F The parameters are further illustrated in 
Figure 6.12. 

E is the energy where the slope of the product curve is zero and 
8 

may be obtained by differentiating equation (6.9) 


( 6 . 10 ) 
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Figure 6.12 The particle spectrum appreciably varies over the threshold region. The effective 

threshold energy, E can deviate from the center. Shown is the power law spectrum, 
efficiency curve and product curve from which the formulation given in the report 
enables Em to be derived. 
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is the effective energy threshold of the channel and is defined 

to be half of the value of the product curve at E for X > — — — , or at 
E s E 2 - E l 

E 2 for X > - — , hence 

E T X ( W = 1/2 (6.11) 

The equation is solved for E Given E^ = 20, E^ = 50, and x = 1 
the value of the product function at E ^ is half the value at E 2 or 

E f ' 1 {E t - 20) = 1/2 30" 1 (30-20) = 1/6 (6.12) 

giving E = 24 keV . Since the center energy, E , equals 25 keV Lhe 
correction is 1 keV at x = 1. For x = 2 the correction is 2.1 keV. 

The spectrum is actually unknown to begin with, requiring that the 
value of x be determined. For initial analysis the threshold, E^, for 
each channel, is taken to be the center energy, E q . The x of the spectrum 
is then derived and used to compute the new E^. Relaxation in this manner 
converges to the proper x in about two itterations. 

Another important consideration is the linearity of the detector to 
different flux intensities. One of the predominant problems is the 
pile-up of output pulses, which simply means the superposition of two or 
more coincident particles giving an erroneous energy count. Each pulse 
is stretched to approximately 5 ysec duration. The probability then that 
a second pulse distorts the first in a random generation in a one second 
interval is 1 divided by 10 + ^ sec = 10"^. Tins occurs because any succeed- 
ing particle within plus or minus 5 ysec appreciably distorts the origi- 
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nal pulse. Allowing three pulses to originate within a one second inter- 
val gives a probability of 2 times 10 ^ for the last two pulses to land 
on the first and 10 for the third to land on the second totaling 3x10 ^ . 
Continuing in this manner for four pulses reveals (3 + 2 + 1)10 ^ = 6x10 ^ 
chances of distorting and so forth. For N pulses 


P = (1 «■ 2 + 3 + . . . + - 1)10 


(6.13) 


where P is the probability that at least one of the // pulses is coinci- 
dent. By equating the opposite terms of the series, the above equation 


reduces to 


atm , 0 ‘ 5 • !L 10 - 5 


(6.14) 


For P = 1/2 the number of particles need only be 320 for at least one 
of the N pulses to be distorted half of the time. 

Figure 6.13 illustrates the detector linearity for channels one 
through four at 70 keV and channels one through three at 20 keV. The 
center energies of the thresholds are 12 keV , 18 keV , 26 keV , and 55 keV 
for channels one through five, respectively. Channel 1 shows the noise 
contribution of 65 ramps per second added to the monitor counts. Also, 
the count rate of channel one is observed to decrease at high rates 
because of the built-in deadtime of 25 ysec. The experimental data 
indicate that the instrument is extremely linear. 

6.3.3 Pulse height analyzer calibration. Two methods were employed 
to obtain pulse height measurements from the accelerator electrons. The 
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Detector Counts x 16 (romps/s) 



Monitor Counts (Counts/s) 


Figure 6.13 Illustrated is the detector linearity curves (experimentally 
derived) for channels one through four at 70 keV and channels 
are through three at 20 keV. Channel 1 shows the noise contri- 
bution of 65 ramps per second added to the monitor counts. 
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first method is the simplest and requires only the measurement of an 
oscilloscope photograph output pulse. The second method used the inte- 
gration of many pulses into pulse height channels. The channel counts 
can be displayed giving the spectrum of the pulses. The peak of the 
spectrum is the desired electron energy. 

The first method is directly accomplished by placing a line 
centrally through the zero line of the noise and another parallel to 
the first through the maximum intensity point of the peak. The difference 
between the lines is the voltage calibration desired for the known electron 
energy. Figure 6.14 is a reproduction of recorded wave shapes for various 
energetic electrons. The final calibration curve relating the voltage 
peak associated with a particular energetic electron is shown in Figure 
6.15. 

The second method requires the use of a pulse height analyzer. A 
sequence of height spectrums are made for particular particle energy 
inputs. Three such spectrums are illustrated in Figure 6.16 for input 
energies of 20, 30, and 40 keV electrons. Again, results can be used to 
determine the final calibration curve for the output voltage peak 
verses incident radiation. 

6.3.4 Sensitivity of the detector to protons. Protons are attenua- 
ted to a large degree because of their inherent mass difference and 
quantum interaction efficiency. The dead zone region consequently reduces 

the proton energy approximately lOkeV before contribution is made to the 

dE 

current pulse. For calculations of the proton range, -r- and depletion 
equivalent energy depth for silicon reference should be made to Figures 
4.7, 4.8 and 4.15, respectively. 
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Figure 0.14 Output pulses from thv preamplifier for various monoenorgetic 


beams of electrons incident on the detector. The height of the 
pulse provides excellent calibration over the entire energy 
range of interest. The upper left hand record shows the noise 
(vertical scale of 50 mV/lioxJ. Directly beneath and adjacent 
is shown the waveforms for 15, 20, 50, 40, 50 and 60 kcV (vertical 
scale of loo mV/ Box J. The next four correspond to energies of 
70, SO, 90 and 100 keV (vertical scale of 200 mV/ Box), finally, 
saturation is shown in the lower right hand signal for a 150 hc\ 
electron (vertical scale of 500 mV/ Box) . 
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Figure 6.15 Final calibration curve for energetic electrons relating the preamplifier 
and postamplifier peak output pulse with the incident particle energy. 
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7. CONCLUSION 

Energetic electrons play a significant role in the upper E-region 
ionization as shown by direct particle measurements and other indirect 
indicators. Analysis of electron-density profiles in the upper E region 
near midnight at Wallops Island indicate that the ionization rate is very 
strongly correlated with geomagnetic activity. No additional source of 
ionization from scattered solar UV radiation is needed to account for 
the ionization rates at Wallops Island. In fact, the present data indi- 
cate that an upper limit for ionization sources other than energetic 
electrons is 10 ^ cm ° sec ^ in the upper E region near midnight. The 
pitch angle distribution of particles and the linear particle count 
rate with altitude suggest that there is a strong concentration of 
particles near and somewhat greater than 90 degrees. This accounts 
for increased ionization in the upper E region give a particular energy 
spectrum since particles are not able to penetrate as low in altitude. 

The spectrums are related by power laws with exponent approximately between 
1 and 4. 

The intermediate layer, observed at night, is explained in terms of 
ionization redistribution caused from horizontal neutral winds. These 
winds induce a converging vertical ion drift which builds up ionization 
at the wind shear (reversal) interface. Using a typical wind pattern 
the continuity equation is solved numerically yielding the modeled layer. 
The layer, using appropriate production, recombination and diffusion 
coefficients, is in good agreement with observed measurements. This 
analysis provides an explanation for why the intermediate layer is more 
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noticeable during quiet geomagnetic conditions than under disturbed con- 
ditions. A wind sensitivity analysis shows that small values of vertical 
ion drift (about 10 m s"*) are sufficient to give rise to a marked inter- 
mediate layer. Applying the continuity equation, given a particular elec- 
tron-density profile, an averaged production rate can be deduced inde- 
pendent of diffusion and of the particular wind system configuration. The 
continuity equation may likewise be solved to determine the winds given 
the intermediate electron-density layer, recombination, and production 
variation. The wind system is observed to descend approximately 1.5 km 
minute'* during the night and can be associated with the tidal winds of 
the (2, 4) -mode of the solar semidiurnal tides. 

The design of a rocket-borne solid-state particle spectrometer was 
undertaken to further enhance the information on the role of energetic 
electrons in the midlatitude ionosphere. The theory, design, calibration 
and operation of this experiment are described in detail. The experiment 
has successfully been flown on four Nike Apache flights, one on 18 April 
1974, and three on 29-30 June 1974, providing plentiful particle data. 

Some of the preliminary results of 18 April 1974 are illustrated in Figures 
2.5 and 2.6. Complete analysis of the above flights will be referred to a 
later report, suffice it to say that flight chart records indicate similar 
behavior to that anticipated in the developed theory of Chapters 2 and 3. 

Future analysis will investigate the ionization relationship in mid- 
latitudes to the anisotropic pitch angle distribution and particle spec- 
trum?. A unified theory needs to be developed relating the satellite 
observations with rocket altitude data and atmospheric interactions. 
Furthermore, the interaction responsible for perturbing the radiation 
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belts and continually replenished then with trapped particles needs 
to be established. Questions relating to geographic variation can be 
assimulated as better interpretations of satellite data are made to 
cohere with an ionization model. This material must be supplemented 
with further rocket investigations to verify the completeness of the 
developed theories. 
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APPENDIX I 


COMPUTER PROGRAM USED TO CALCULATE IONIZATION RATES 
FROM PARTICLE FLUX AND SPECTRUM INFORMATION 
The program used to compute electron ionization altitude profiles 
of Figure 2.12 was based upon the methods enunicated by Wulff and 
Gledhill [1973]. The pitch angle distribution, employed in this develop- 
ment, is assumed to be isotropic over the downward hemisphere. Although 
this approximation is not necessarily valid the resulting attenuation 
curves are extremely helpful in determining overall sensitivity of 
profiles to spectrum and interaction coefficients. This program was, in 
the main, developed by Monro [1974]. The height profiles are 
calculated for exponential spectrums with exponents of 1.0 thru 7 in 
half interval steps. The relevant ionization of each major constituent 
(0 + , C>2 + , N^ + , N + ) is printed verses the altitude. The input data from 
Table 1 of Wulff and Gledhill [1973] is read into array A(20) . The 
number densities are read into AN (4, 20). 

c production ratsS f^om energetic electrons 

DIH£»SJoN A (20) , k T < 2«* ) »AN< 4,20) .2(20 ).Q(5) 

cona s lugf (10 . ) 

HE'D i»n»m 

1 T * ? I 5 ) 

READ ?» <M I >• l*l»M) 

2 F ORmAT(iOF0.O ) 

READ 2*(MTCI).I«1*N) 

DO 3 I *1*4,1 

READ 2 * C *N < I * J) , J*1»M 

DO 3 J S 1 » N, 1 

3 AN<I*J)*fcXPF*CAN< I,J)*CONA) 

READ 4 • < Z( I ) » I *1, N) 

4 FORMAT l gElO . 0 ) 

READ 2* ECh I N»EOMA x » DH) 

E0=E0hIn 

13 PRINT I 5 .FC 

15 FORMAT(2*,4HfcO * * F5 , 1 / ) 

PRinI 1& 
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16 ro^M^T t 4X # 6MHg l G* T, 7* , 3*00* » 8X » 4HQ02 4 * 8X* 

14H0N2* »9X,3HQN*«8X*4HUT0l) 

DO 12 1*1 ,N* 1 
0(5 ) *0 • 

6N*1*0 

10 RP*4*17e- 6«EN**1.74 
ZRP*Z( 1 )/PP*l. 064 
IF< ZPP'l « 2 >5* 5 * 6 

5 J*1 
P*. t 

9 IF < ZPP”P ) 7» 7# 8 

8 JiJ+1 

P*P*. 1 
GO TO 9 

7 AFT 8 M j* 1 )4<P-ZR5)/.1*(A( J)-A( J*l) ) 

AZRsA r -T # EN*l. 63/^9 
R*En/EO 

F4saZR*A^< 4, I )»10 0./PXPF(R)/35./60/l.E3 
Q(5) = Qt5>*l : H 

6 E N s c 'W . i 

I F t -2 q . >10*10*11 

11 I r <Q<5>-l.E-12) 23 .20 **21 

20 DO ?/c Jsl , 5 * 1 

22 Q C J ) 5 0 • 

GO TO I? 

21 B< S 1**1 # 15«(AN(2» 1 > ♦ 4 N ( 1 # I )/2t)/AM(3. I ) 

0(l) s ( . 2 9 *aN ( 2 * 1 ) ♦AM < 1 . I )/2. )»1.15/AN(3. I ) *0 ( 5 ) / 9K 
0<2>*i.t5#.71*AN<2* I l/AM(3. I >/PK«G(5> 
0(3js.bt*(0(5)-0tl > - Q < 2 > » 

Q(<) = .!9«(0(5)-gm-G(?n 
1? P R I m T 17 »HT ( I ) « (Q( J) • J = l*5) 

17 FDh m at(2X,F7.2,2X.5E12.3) 

E3 s F>-'*DgO 

I?(gL.EoMAX)18»le*lA 
IP PR ( mT 1 q 
19 FDH M ATliHo////) 

GO TO 13 
14 CAUL E* t T 
END 


REPRODUCIBILITY of the 

ORIGINAL page is poor 


J 



1 


! 


198 


APPENDIX II 

SOLUTION OF TIME -INDEPENDENT CONTINUITY EQUATION WITHOUT DIFFUSION 

Two methods of solution were persued; 1) an e.nalytical representation 
and 2) a numerical computer algorithm. The appropriate continuity equation 
to be solved is 

N ^* w l^ =q ~ aN2 (A2.ll 

where w is the vertical ion wind taken to be sinusoidal (equation 3.8) 
and q and a are the production and recombination coefficients assumed to 
be constant for the analytical representation. The equation, in general, 
is nonlinear with variable coefficients. 

Analytical solution. Equation (A2.1) is classified as a form 
of the Riccati equation given as 

+ Qiz)y * RWy 1 2 • P{ 3) CA2.2) 

For the special case of q = 0 the substitution of y • ^ transforms the 
equation into a first order linear equation which can be directly solved. 

The nonlinear first order equation (A2.2) can be transformed into a 
second order linear differential equation 

2 

R ~ - (/?' - gy - B 0 (A2.3) 

dz z 

where y is given by 


1 du 

V " R: 7 3F 


(A2.4) 
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For the case of interest a and q are constant and dw/dz is given by 


3 3 ^ _ . z -v 2 tt 

w = ~~ (20 sm— tt) = — tt cos^r z 

3 z 2 2 v 30 3 30 


(A2.S) 


proceeds as follows. The appropriate R s Q> and P are 


D a ot z 

R ~ “7 = TR csc t77 77 

W 20 30 


^ 1 dh) 7T Z 

Q - — -z — = Yn ctn t TT 77 

U 3s 30 30 


(A2.6) 


P = 2- = -S. CSC JL tt 
w 20 30 


Making these substitutions into equation (A2.3) and rearranging gives 


° 2 + (Ctn if " + Ctn if ^ 0 


q a 


400 sin 2 tt 


M = 0 


where D is the differential operator. 

C0 STT2 

Now letting X = — , the reduced equation is formulated as 


D 


2 


+ aD 


b 

X{l-X 2 ) 


u = 0 


(A2 .7) 


where 

2 

IT J r <?a 

“ ■ 450 and b ■ 400 

Equation (A2.7) is a linear second order differential equation. The 
state space A matrix is given as 


X = 


0 1 

b 

2 a 
X(1 - X ) 


(A2.8) 
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which is applicable to standard analysis methods. Also equation (A2.7) 
may be solved using the WKB method of solution. 

Numerical Solution. The equation (A2.1) was found to be 
quite sensitive requiring every precaution to insure stability. A 
predictor-corrector method of the Euler-Romberg type was found to operate 
satisfactorily. Further reference to the theory of this method is found 
in McCalla [1967], The program used to solve the continuity equation, 
employing this technique, is provided on the following pages. 


i 



i , 


T 


’"'T f' 


T" 


201 


FIR TP AN PROGRAM FOR clJLEP-8CM° ERG METHOD 

■HI.:iL^SI CEDU1QQU Yi.lQQlt T.l 2QJ 


CrMMCN LIMIT 
^ 1 = 1 , 1 ^ 15 S 21 L 


C.. 


X. 

C 


L I M IT=1 

_ IN I T LR- V A L'JIS . . STEP _.SI 7 C»...CQN'VFRGENCE. TERM,, 
P540 (5.100 X ( 1 ) • Y ( 1 ) * H, EPS. L M 4X » NC 
CQ Mpi)T = niju = p_j£ AL SH L UTI ONJ/Al UF. FO R EACH STE° _N 


COMMUTE VALUES »Y EULER MfTMCO 

. A . 14*1 

7 GO Tr. « 

._e i._li. m lt =: 

9 X(N+i)=X(N)*w 

„J_Q. Xi’ifcUs.UU. 



N=M+r 

5 _T.i 11- V t N1 +H* PC T ( X ( N ) , Y ( N) ) 


IF (LIMIT) 
_&_L = 1 


1 3 » 5 


10 XC=X(N) 
jm=y t.Ni. 


17 


19 

-20 


L I M=?**L 

.00 .15. 

F = FCT (XC.YC) 

_L£__U_ IM1T ) -2A*2_. 19 


2 i 
? 2 
23 


14 XC = XC + H/Fl.OATU imj 

1 5 YC=YC-HH/Ft.nAT ( L 111) 1*£ 

T(L+1)=YC 

C COMPUTE. £LLL£R_F-C M 9 IF. G_XA 3.L E- VALUES-UKl 

’4 M = 1 

— 7-5 JU. L - 


26 


28 
-O i. 


20 T(K )=(2.*»m«t(K+1)-T(k))/( 

IF.. 1.K.-L1, . 2 5 1 35.. 25 

25 I c (ARS(T(K)-T(K-l) )-EPS) 50,50.30 
-40- -W—M*! 


K = K— 1 


32 

35 

IF (L-LMAX) 40,75. 

1 “3 

4H 

i =1 4-1 

3 H 


GO TO 10 

. _^_5 

40- 

XI M*1 ) =X ( w ) +H 

36 

55 

Y ( M ♦ 1 ) = T ( K ) 

. — 17 


S=FrT(X(N).Y(N)) 


-j-a- 


XS-X(N)»l.E-5+120. 

■UsS ILK )»20. 


40 WRITE (6,110) X£,Y(N),S, W 

-.1 tR= 4 U 4 

42 IF(N-MC) 5 , oC , 6 0 

— L tUimiL SnillTIPN P 0 T M T R f.nMPII TFn AY EU LER R QM B£R G 

93 60 WRITE (6.105) 

-44. gn AS „Nsl .NC 


/ ^ 


S = FCT ( X ( N ) ,Y(M) ) 


L 20 . REW^pUC^ILjTY pFjJiE 
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-J 


M 


t 

..J 


<^7 * = $ 1*4(01/3 ,E4-6»*X(N) )*20. 

— 4 J k&XI £ . .ItnUai XS > Y t NJuSo-Jtl — 

49 continue 

— SO. on T n jja 

r ptjTO'JT SP.UJTIfl* POINTS PP I QR TO CONVERGENCE FAILURE 

— SX 1 2S_4?a 7S_liu»i lit 

5 2 N?«N 

WiUXP -16. 110^ IX-IN)* Y(Ni -N*l«N2J 

50 STOP 

I fl fl -F.-3 F 2& L Ljp.1^7^ : FlflxJU. LI II — — 

105 e-DR VAT ( 2X. IGMT/KM , 8x. 1 2H0ENS ITY/C M -3 *3X 1 1 5H0ER I VI TI VF/C M -« ♦ » ?X 

tLLMhlKQ/ M/SEC ) 

no F0P M AT ( IX * FlO.f.TXtElZ.AtSXtFlO.Sf^XfFlO.^) 
ilO.fOlLVAI. I20H .FAILURE ID CONVERGE) 

tMC 


FUnCT I qN FCT U RGI,ARG2) 
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_ .3 FCl*fl . 

RETURN 
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APPENDIX III 

COMPUTER PROGRAM USED TO SOLVE THE TIME- INDEPENDENT 
CONTINUITY EQUATION WITH DIFFUSION 

The steady-state continuity equation can be written as 


_d 

dz 


(- 0 >/.’) = q - 


a N 2 + 


sin I 


d 

3F 




(A. 0 


where z is the altitude, w is the ion vertical wind velocity, is the 
electron density, q is the ionization rate, a is the recombination coeffi- 
cient, I is the dip angle and D is the plasma (or arbipolar) diffusion 
coefficient. The particular wind model employed for the converging ion 
wind is a sinusoid 


w 


20 sin 


TT (150 - z) 
30 


(A3. 2) 


For initial analysis D, q, and a were taken to be realistic constants 
appropriate to 150 kilometers. 

The equation is very sensitive and prone towards instability if 
slight errors are encountered. Difference equation methods for solving 
the equation resulted in ir.stab : ity in every case tried, even for 
implanted damping factors. 

An iterative integration method was employed which gave promising 
solutions; however, required much computer time for analysis. This 
method seems to be applicable to a much more generalized analysis and 

therefore outlined. Integration of equation (A3.1) from 0 to z yields 

dN 

a first order differential equation where u = ^ = 0 at 2 = 0, namely 


i 


t 




dN w„ 1 f * .,2, . , A - 

Is ' D N D \ (<? - ( A3 - 3 ) 

J Q 

Now, any first order equation has a solution of the form 

f 3 

P(s)d3 

J f" 

N = e ° Q(z) e 

* o 

W l Z 2 

where P(s) is and $(s) is — ((7 - aiV )dz . 

The procedure is to assume an /V distribution and thereby determine 
4 ( 2 ) and then N[z) by application of (A3. 4). This value of N{z) is again 
used to calculate a new Q(z ) and N(z ) . Relaxation converges upon the 
desired profile if sufficient accuracy is taken. 

The best method found for solving equation (A3.1) is to employ the 
fourth order Runge-Kutta method using 0.5 km steps. The procedure requires 
specification of an initial slope and magnitude as a starting boundary con- 
dition, however, the actual boundary conditions known are zero slope at 
the peak and valleys of the profile. To make these boundary conditions 
compatible with those demanded by the Runge-Kutta method a relaxation 
process was employed. The initial slope was set at zero and a zero 

diffusion value was initially picked as the starting value. The marching 
solution then diverged either positively or negatively. If positive the 
initial value was reduced and vice versa until the condition of zero 
slope was rendered at the peak and valley. Usually five to twelve iterations 
are required for the final solution to lie within the chosen limits. This 
procedure is reproduced on the following pages. 


•2 

P{z)dz 

'o dz (A3. 4) 
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0060 
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12 YU= Y ( 1 J 
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GO tq s 

. ...0064. 
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0065 
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0066 
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21 i=n-i _ . . 
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APPENDIX IV 


COMPUTER PROGRAM USED TO DERIVE WINDS FROM AN ELECTRON-DENSITY PROFILE 
To calculate the winds, w(z) , from an electron-density profile, 

N{z), the derived equation (3.9) applies 


w(3) = mr {z ■ z \ )q ' aHz) & 


(A4.1) 


where the appropriate q is used for the lower and upper portion of the 
profile. The average ionization q is first calculated in the program by 


use of Simpson's rule. 


Z 2 ’ 3 1 


‘ 2 JIT 2 J 

aN az 


The recombination coefficient a is taken to be variable as discussed 


in Section 3.2.2. 
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